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The nonlinear rheological properties of dense colloidal suspensions under steady shear are 
discussed within a first principles approach. It starts from the Smoluchowski equation of in- 
teracting Brownian particles in a given shear flow, derives generalized Green-Kubo relations, 
which contain the transients dynamics formally exactly, and closes the equations using mode 
© , coupling approximations. Shear thinning of colloidal fluids and dynamical yielding of col- 

. loidal glasses arise from a competition between a slowing down of structural relaxation, be- 

cause of particle interactions, and enhanced decorrelation of fluctuations, caused by the shear 

Sadvection of density fluctuations. The integration through transients approach takes account of 
the dynamic competition, translational invariance enters the concept of wavevector advection, 
and the mode coupling approximation enables to quantitatively explore the shear-induced sup- 
pression of particle caging and the resulting speed-up of the structural relaxation. Extended 
q comparisons with shear stress data in the linear response and in the nonlinear regime measured 

i in model thermo-sensitive core-shell latices are discussed. Additionally, the single particle mo- 

tion under shear observed by confocal microscopy and in computer simulations is reviewed 
and analysed theoretically. 
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List of Abbreviations and Symbols 

Goo Shear modulus of a solid (transverse elastic constant or Lame-coefficient) 
rjo Newtonian viscosity of a fluid 
a Shear stress 
7 Shear rate 

g(t) Time dependent shear modulus; in the linear response regime denoted as g lr (t) 
of the quiescent system; generalized one if including dependence on shear rate 
r Maxwell (final or a- relaxation) time of structural relaxation 
G' (u>) Storage modulus in linear response 
G"(u) Loss modulus in linear response 
r/ Shear viscosity; defined via r\ — cr(j)/j 
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s q 


Equilibrium structure factor 


l_) 

Hh 


Hydrodynamic radius of a colloidal particle (radius a = Rh assumed) 


a 


Colloid diameter (d = 2a = 2Rh assumed throughout) 


Kb -i 


1 hermal energy 




Solvent viscosity 


J->0 


Stokes Einstein Sutherland diffusion coefficient at infinite dilution 


Peo 


Bare Peclet number 


Pe 


Dressed Peclet or Weissenberg number 


<P 


Packing fraction <j> = ^-R^ H n of spheres of radius Rh at number density n 


e 


Separation parameter in MCT giving the relative distance in a thermodynamic 




control parameter to its value at the glass transition 


A 


MCT exponent parameter 


Goo 


Instantaneous isothermal shear modulus 


?7oo 


High frequency viscosity 


a+ 


Dynamic yield stress of a shear molten glass 


HI 


Hydrodynamic/ solvent induced interactions 


MCT 


Mode coupling theory 



ITT Integrations through transients approach 
SO Smoluchowski operator fl 

PY Percus-Yevick theory giving the approximate PY S q of a hard sphere fluid 
ISHSM Isotropically sheared hard spheres model 

F{2 Schematic model without spatial resolution considering a single correlator 
1 

Introduction 

Rheological and elastic properties under flow and deformations are highly charac- 
teristic for many soft materials like complex fluids, pastes, sands and gels, viz. soft 
(often metastable) solids of dissolved macromolecular constituents [1]. Shear de- 
formations, which conserve volume but stretch material elements, often provide the 
simplest experimental route to investigate the materials. Moreover, solids and fluids 
respond in a characteristically different way to shear, the former elastically, the latter 
by flow. The former are characterized by a shear modulus Goo, corresponding to a 
Hookian spring constant, the latter by a Newtonian viscosity 770, which quantifies the 
dissipation. 

Viscoelastic materials exhibit both, elastic and dissipative, phenomena depend- 
ing on external control parameters like temperature and/ or density, and depending on 
frequency or time-scale of experimental observation. Viscoelastic fluids differ from 
pastes and sands in the importance of thermal fluctuations causing Brownian motion, 
which enables them to explore their phase space without external drive like shak- 
ing, that would be required to fiuidize granular systems. The change between fluid 
and solid like behavior in viscoelastic materials can have diverse origins, including 
phase transitions of various kinds, like freezing and micro-phase separation, and/or 
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molecular mechanisms like entanglement formation in polymer melts. One mecha- 
nism existent quite universally in dense particulate systems is the glass transition, 
that structural rearrangements of particles become progressively slower [2] because 
of interactions/ collisions, and that the structural relaxation time grows dramatically. 

Maxwell was the first to describe the viscoelastic response at the fluid-to-glass 
transition phenomenologically. He introduced a time-dependent shear modulus g(t) 
describing the response of a viscoelastic fluid to a time-dependent shear deformation, 



Here a is the (transverse) stress, the thermodynamic average of an off-diagonal 
element of the microscopic stress tensor, and j(t) is the time-dependent shear 
rate impressed on the material starting at time t = 0. Maxwell chose the Ansatz 
g(t) = Goo ex P {~(t/ T )}> which interpolates inbetween elastic behavior a(t — > 
0) ~ Goo j(t) for short times t <C r and dissipative behavior, a(t) w r/ j(t), for 
long times, i»r; the strain 7(i) is obtained from integrating up the strain rate, 
j(t) = f*dt'j(t'). Maxwell found the relation 770 = G m t which connects the 
structural relaxation time and the glass modulus Goo to the Newtonian viscosity. He 
thus explained the increase of the viscosity at the glass transition by the slowing 
down of the structural dynamics (viz. the increase of t), and provided a definition of 
an idealized glass state, where r = 00. It responds purely elastically. 

Above relation (1) between a and 7 is exact in linear response, where non-linear 
contributions in 7 are neglected in the stress. The linear response modulus (to be 
denoted as g h (t)) itself is defined in the quiescent system and describes the small 
shear-stress fluctuations always present in thermal equilibrium [1,3]. Often, oscil- 
latory deformations at fixed frequency lu are applied and the frequency dependent 
storage- (G'(w)) and loss- {G"(io)) shear moduli are measured in or out of phase, 
respectively. The former captures elastic while the latter captures dissipative contri- 
butions. Both moduli result from Fourier-transformations of the linear response shear 
modulus g h (t), and are thus connected via Kramers-Kronig relations. 

The stationary, nonlinear rheological behavior under steady shearing provides 
additional insight into the physics of dense colloidal dispersions [1,3]. Here, the 
shear rate is constant, j(t) = 7, and the stress in the stationary state achieved after 
waiting sufficiently long (taking t — > 00 in Eq. (1)) is of interest. Equation (1) may 
be interpreted under flow to state that the non-linearity in the stress versus shear 
rate curve (the relation 17(7) is termed 'flow curve') results from the dependence of 
the (generalized) time-dependent shear modulus g(t,j) on shear rate. The (often) 
very strong decrease of the viscosity, defined via 77(7) = a(j)/j, with increasing 
flow rate is called 'shear thinning', and indicates that the particle system is strongly 
affected by the solvent flow. One may thus wonder whether the particles' non-affine, 
random motion relative to the solvent differs qualitatively from the Brownian motion 
in the quiescent solution. Taylor showed that this is the case for dilute solutions. A 
single colloidal particle moves super-diffusively at long times along the direction 
of the flow. Its mean squared non-affine displacement grows with the third power 
of time, much faster than the linear in time growth familiar from diffusion in the 
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quiescent system . A priori it is thus not clear, whether the mechanisms relevant 
during glass formation in the quiescent system also dominate the nonlinear rheology. 
Solvent mediated interactions (hydrodynamic interactions, HI), which do not affect 
the equilibrium phase diagram, may become crucially important. Also, shear may 
cause ordering or layering of the particles leading to heterogeneities of various kinds 
[4]. 

Within a number of theoretical approaches a connection between steady state 
rheology and the glass transition has been suggested. Brady worked out a scaling de- 
scription of the rheology based on the concept that the structural relaxation arrests at 
random close packing [5]. In the soft glassy rheology model, the trap model of glassy 
relaxation by Bouchaud was generalized by Cates and Sollich and coworkers to de- 
scribe mechanical deformations and ageing [6-8]. The mean field approach to spin 
glasses was generalized to systems with broken detailed balance in order to model 
flow curves of glasses under shear [9, 10]. The application of these novel approaches 
to colloidal dispersions has lead to numerous insights, but has been hindered by the 
use of unknown parameters in the approaches. 

Dispersions consisting of colloidal, slightly polydisperse (near) hard spheres ar- 
guably constitute one of the most simple viscoelastic systems, where a glass tran- 
sition has been identified. It has been studied in detail by dynamic light scattering 
measurements [11-19], confocal microscopy [20], linear [21,22], and non-linear 
rheology [23-32]. Computer simulations are available also [33-35]. Mode coupling 
theory (MCT) has provided a semi-quantitative explanation of the observed glass 
transition phenomena, albeit neglecting ageing effects [36] and decay processes at 
ultra-long times that may cause (any) colloidal glass to flow ultimately [2, 37, 38]. It 
has thus provided a microscopic approach recovering Maxwell's phenomenological 
picture of the glass transition; and r could be calculated starting from the par- 
ticle interactions as functions of the thermodynamic control parameters. MCT was 
also generalized to include effects of shear on time dependent fluctuations [39-41], 
and, within the integrations through transients (ITT) approach, to quantitatively de- 
scribe all aspects of stationary states under steady shearing [42-44]. 

The MCT-ITT approach thus provides a microscopic route to calculate the gener- 
alized shear modulus g(t, 7) and other quantities characteristic of the quiescent and 
the stationary state under shear flow. While MCT has been reviewed thoroughly, see 
e.g. [2, 37, 38], the MCT-ITT approach shall be reviewed here, including its recent 
tests by experiments in model colloidal dispersions and by computer simulations. 
The recent developments of microscopy techniques to study the motion of individual 
particles under flow and the improvements in rheometry and preparation of model 
systems, provide detailed information to scrutinize the theoretical description, and to 
discover the molecular origins of viscoelasticity in dense colloidal dispersions even 
far away from thermal equilibrium. 

1 This effect that flow speeds up the irreversible mixing is one mechanism active when stir- 
ring a solution. The non-affine motion even in laminar flow prevents that stirring backwards 
would reverse the motion of the dissolved constituents. 
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The outline of the review is as follows: At first, the microscopic starting points, 
the formally exact manipulations, and the central approximations of MCT-ITT are 
described in detail. Section 3 summarizes the predictions for the viscoelasticity in 
the linear response regime and their recent experimental tests. These tests are the 
quantitatively most stringent ones, because the theory can be evaluated without tech- 
nical approximations in the linear limit; important parameters are introduced here, 
also. Section 4 is central to the review, as it discusses the universal scenario of a 
glass transition under shear. The shear melting of colloidal glasses and the key phys- 
ical mechanisms behind the structural relaxation in flow are described. Section 5 
builds on the insights in the universal aspects and formulates successively simpler 
models which are amenable to complete quantitative analysis. In the next Section, 
those models are compared to experimental data on the microscopic particle motion 
obtained by confocal microscopy, to data on the macroscopic stresses in dispersions 
of novel model core-shell particles close to equilibrium and under steady flow, and 
to simulations providing the same information for binary supercooled mixtures. In 
the last Section, recent generalizations and open questions are addressed. 

2 

Microscopic approach 

MCT considers interacting Brownian particles, predicts a purely kinetic glass tran- 
sition and describes it using only equilibrium structural input, namely the equilib- 
rium structure factor S q [3, 45] measuring thermal density fluctuations. MCT-ITT 
extends this Statistical Mechanics, particle based many-body approach to dispersions 
in steady flow assuming a linear solvent velocity profile, but neglecting the solvent 
otherwise. 



Interacting Brownian particles 

N spherical particles with radius Rh are considered, which are dispersed in a vol- 
ume V of solvent (viscosity r/g). Homogeneous shear is imposed corresponding to 
a constant linear solvent velocity profile. The flow velocity points along the x-axis 
and its gradient along the y-axis. The motion of the particles (with positions (t) for 
i = 1, . . . , N) is described by N coupled Langevin equations [45] 



Solvent friction is measured by the Stokes friction coefficient £ = GirrjsRn- The 
interparticle forces Fj = —djdvi U({rj}) derive from potential interactions of par- 
ticle i with all other colloidal particles; U is the total potential energy. The solvent 
shear-flow is given by v solv (r) = 7 y x , and the Gaussian white noise force satisfies 
(with a, (i denoting directions) 



2.1 




(2) 
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(f?(t) //(*')> - 2C fc B T5 a/3 5(i - t') 



where ksT is the thermal energy. Each particle experiences interparticle forces, sol- 
vent friction, and random kicks from the solvent. Interaction and friction forces on 
each particle balance on average, so that the particles are at rest in the solvent on 
average; giving for their affine motion: (r, (i)) = r j (0) + 7 1 yt (0) x . The Stokesian 
friction is proportional to the particle's motion relative to the solvent flow at its posi- 
tion; the latter varies linearly along the y-direction. The random force on the level of 
each particle satisfies the fluctuation dissipation relation. The interaction forces Fj 
need 

Even though Eq. (2) thus has been obtained under the assumption, that solvent 
fluctuations are close to equilibrium, the Brownian particle system described by it 
may reach macroscopic states far from thermal equilibrium. Moreover, under (finite) 
shear, this holds generally because the friction force from the solvent in Eq. (2) can 
not be derived from a conservative force field. It has non-vanishing curl, and thus the 
stationary distribution function \P describing the probability of the particle positions 
r-j can not be of Boltzmann-Gibbs type [46]. 

Already the microscopic starting equation (2) of MCT-ITT carries two important 
approximations. The first is the neglect of hydrodynamic interactions (HI), which 
would arise from the proper treatment of the solvent flow around moving parti- 
cles [3,45]. As vitrification is observed in molecular systems without HI, MCT-ITT 
assumes that HI are not central to the glass formation of colloidal dispersions. The 
interparticle forces are assumed to dominate and to hinder and/or prevent structural 
rearrangements close to arrest into an amorphous, metastable solid. MCT-ITT as- 
sumes that pushing the solvent away only provides some additional instantaneous 
friction, and thus lets short-time transport properties (like single and collective short- 
time diffusion coefficients, high frequency viscoelastic response, etc.) depend on HI. 
The second important approximation in Eq. (2) is the assumption of an homogeneous 
shear rate 7. This assumption may be considered as a first step, before heterogeneities 
and confinement effects are taken into account. The interesting phenomena of shear 
localization and shear banding and shear driven clustering [47-51] therefore are not 
addressed. All difficulties in Eq. (2) thus are connected to the many-body interactions 
given by the forces Fj, which couple the N Langevin equations. In the absence of 
interactions, F, = 0, Eq. (2) immediately leads to the super-diffusive particle motion 
mentioned in the introduction, which often is termed 'Taylor dispersion' [45]. 

As is well known, the considered microscopic Langevin equations, are equivalent 
to the reformulation of Eq. (2) as Smoluchowski equation; it is a variant of a Fokker- 
Planck equation [46]. It describes the temporal evolution of the distribution function 
^{{vi] , t) of the particle positions 



d t *{{v i },t) = f2*{{T i },t) 



(3a) 



employing the Smoluchowski operator (SO) [3,45], 




(3b) 



8 



Matthias Fuchs 



built with the Stokes-Einstein-Sutherland diffusion coefficient D = k B T/( of a 
single particle. Averages performed with the distribution function <F agree with the 
ones obtained from the explicit Lagevin equations. 

The Smoluchowski equation is a conservation law for the probability distribution 
in coordinate space, 



formed with probability current j. Stationary distributions, which clearly obey d t & s — 
0, which are not of equilibrium type, are characterised by a non-vanishing probability 
flux j| 7^ 0, where 



Under shear, j s can not vanish, as this would require the gradient term to balance 
the term proportional to 7 which, however, has a non-vanishing curl; the 'potential 
conditions' for an equilibrium stationary state are violated under shear [46]. 

The ITT approach formally exactly solves the Smoluchowski equation, follow- 
ing the transients dynamics into the stationary state. In this way the kinetic com- 
petition between Brownian motion and shearing, which arises from the stationary 
flux, is taken into account in the stationary distribution function. To explicitly, but 
approximatively compute it, using ideas based on MCT, MCT-ITT approximates the 
obtained averages by following the transient structural changes encoded in the tran- 
sient density correlator. 



Generalized Green-Kubo relations 

Formally, the H-theorem valid for general Fokker-Planck equations states that the 
solution of Eq. (3) becomes unique at long times [46]. Yet, because colloidal particles 
have a non-penetrable core and exhibit excluded volume interactions, corresponding 
to regions where the potential is infinite, and the proof of the H-theorem requires 
fluctuations to overcome all barriers, the formal H-theorem may not hold for non- 
dilute colloidal dispersions. Nevertheless, we assume that the system relaxes into a 
unique stationary state at long times, so that <F(t — > 00) = <F S holds. This assumption 
is self-consistent, because later on MCT-ITT finds that under shear all systems are 
'ergodic' and relax into the stationary state. In cases where phase space decomposes 
into disjoint pockets ('nonmixing dynamics'), the distribution function calculated in 
Eq. (4) averages over all compartments, and can thus not be used. 




_d_ _J_ 



2.2 

Integration through transients (ITT) approach 
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As alread stated, homogeneous, amorphous systems are assumed so that the sta- 
tionary distribution function <F S is translationally invariant but anisotropic. The for- 
mal solution of the Smoluchowski equation for the time-dependent distribution func- 
tion 

>p(t) = e nt <P e (4a) 

can, by taking a derivative and integrating it up to t = oo, be brought into the form 
[42,43] 

*. = *« + r^ / dt# e a xy e nU , (4b) 

where the adjoint Smoluchowski J?' operator arises from partial integrations over 
the particle positions (anticipating that averages built later on with are done by 
integrating out the particle positions). It acts on the quantities to be averaged with 
S's . The assumption of spatial homogeneity rules out the considerations of thermody- 
namic states where the equilibrium system would e.g. be crystalline. The equilibrium 
state is described by >P e , which denotes the equilibrum canonical distribution func- 
tion, )P e oc e~ l 1 1 ( fcsT ), which is the time-independent solution of Eq. (3a) for 7 = 0; 
in Eq. (4b), it gives the initial distribution at the start of shearing (at t = 0). The 
potential part of the stress tensor <r xy = — ^2^ =1 Ff yi entered via f}ip e = 7 <r xy \l/ e . 
The simple, exact result Eq. (4b) is central to the ITT approach as it connects steady 
state properties to time integrals formed with the shear-dependent dynamics. Ad- 
vantageously, the problem to perform steady state averages, denoted by (. . has 
been simplified to performing equilibrium averages, which will be denoted as (. . .) 
in the following, and contain the familiar <P e . The transient dynamics integrated up in 
the second term of Eq. (4b) contains slow intrinsic particle motion, whose handling is 
central to the MCT-ITT approach. Generalized Green-Kubo relations, formally valid 
for arbitrary 7, can be derived from Eq. (4b). 

The adjoint Smoluchowski operator was obtained using in the partial integrations 
over the particle positions the incompressibility condition, Tracc{«} = 0, which 
should always holds for the solvents of interest in this review. It takes the explicit 
form (where boundary contributions are neglected throughout, simplifying the partial 
integrations): 

f2^ =Y / (d l +F l +r l -K T )-d t . 

i 

This formula already uses a handy notation 2 employing the shear rate tensor k = 
7 xy (that is, n a p = jS ax S/3 y ), and dimensionless quantities. They are introduced 
by using the particle diameter d as unit of length (throughout we convert d = 2Rh), 
the combination d 2 /Do as un it of time, and ksT as unit of energy, whereupon the 
shear rate turns into the bare Peclet number Pe = ■jd 2 /D . It measures the effect of 

2 The simplified notation with dimensionless quantities is used in the Sections containing 
formal mainpulations, and in a number of original publications. 
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affine motion with the shear flow compared to the time it takes a single Brownian par- 
ticle to diffuse its diameter d. One of the central questions of the nonlinear rheology 
of dense dispersions concerns the origin of very strong shear-dependences in the vis- 
coelasticity already for (vanishingly) small bare Pe numbers. Thus we will simplify 
by assuming Pe <C 1, and search for another dimensionless number characterizing 
the effect of shear. 

The formally exact general result for \P S in Eq. (4b) can be applied to compute 
the thermodynamic transverse stress, (7(7) = (cr xy ) /V. Equation (4b) leads to an 
exact non-linear Green-Kubo relation: 



a( 7 )= 7 / dtg(t,i), (5a) 







where the generalized shear modulus g(t, 7) depends on shear rate via the Smolu- 
chowski operator from Eq. (3b) 



This relation is nonlinear in shear rate, because 7 appears in the time evolution oper- 
ator J?t, the adjoint of Eq. (3b). In MCT-ITT, the slow stress fluctuations in g(t, 7) 
will be approximated by following the slow structural rearrangements, encoded in 
the transient density correlators. 

But, before discussing approximations, it's worthwhile to point out that formally 
exact explicit expressions for arbitrary steady-state averages can be obtained from 
Eq. (4b). Using the definition /( 7 ) = (/ q=0 )^ /V, where q is the wavevector, and 
/ q=0 = J rfr/(r) denotes the integral over an arbitrary density /(r), one finds the 
general generalized Green-Kubo relation: 

/•OO 

/( 7 ) = (/ q=0 )/y+^ jdt (<J xy e nU Af^o) , (5c) 

where the symbol AX for the fluctuation in X was introduced, AX = X - (X), 
because all mean values (which are constants, for these purposes) drop out of the ITT 
integrals, leaving only the fluctuating parts to contribute. Generalizations of Eq. (5c) 
valid for structure functions (see e.g. Eq. (6b)) and stationary correlation functions 
(see Eq. (8)) are presented in Ref. [43]. Note that all the averages, denoted (...), 
are evaluated within the (Boltzmann) equilibrium distribution <P e . Why only q = 
appears in Eq. (5c) is discussed in Sect. 2.2.2 . 

It is these generalized Green-Kubo relations Eq. (5c) which are formally exact 
even for arbitrary strong flows, and which form the basis for approximations in the 
MCT-ITT approach. These approximations are guided by the evident aspect that slow 
dynamics strongly affects the time integral in Eq. (5c). Therefore, in MCT-ITT ap- 
proximations are employed that aim at capturing the slow structural dynamics close 
to a glass transition. It would be interesting to employ the generalized Green-Kubo 
relations also in other contexts, where e.g. entanglements lead to slow dynamics in 
polymer melts. 
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2.2.2 

Aspects of translational invariance 

The generalized Green-Kubo relations contain quantities integrated/ averaged over 
the whole sample volume. Thus, the aspect of translational invariance/ homogeneity 
does not become an issue in Eq. (5) yet. A system is translational invariant, if the 
correlation between two points r and r' depends on the distance r r' between 
the two points only. The correlation must not change if both points are shifted by 
the same amount. (Additionally, any quantity depending on one space point r only, 
must be constant.) A system would be isotropic, if additionally, the correlation only 
depended on the length of the distance vector, |r — r'|; but this obviously can not be 
expected, because shear flow breaks rotational symmetry of the SO in Eq. (3b). Shear 
flow also breaks translational symmetry in the SO of Eq. (3b), therefore it is a priori 
surprising, that translational invariance holds under shear. Moreover, discussion of 
translational invariance introduces the concept of an advected wavevector, which 
will become important later on. 

The time-dependent distribution function \P(t) from Eq. (4a) can be used to show 
that a translationally invariant equilibrium distribution function <P e leads to a transla- 
tionally invariant steady state distribution \P S , even though the SO in Eq. (3b) is not 
translationally invariant itself. To show this, a point in coordinate space (ri, . . . , rjv) 
shall be denoted by r, and shall be shifted, r — > r', with r'j = r, + a for all i; a is 
an arbitrary constant vector. This gives 

flt(r) = flt(r') - a k t P , withP = ^0 4 , 

i 

explicitly stating that the SO is not translationally invariant. From Eq. (4a) follows 

where \P e (r') = &e(r) was used. The SO Q and the operator P«a commute, 
because the shear rate tensor satisfies k ■ k = 0, and because the sum of all internal 
forces vanishes due to Newton's third law: 

(P«a) Q-Q (P«a) = 
Yl { 9 ^ d i ■ ~ ( 9 i ' fr^ 1 Ka ~ [(^di) {rjK T dj) - (d j Kr j )(d l na)] 



>'j 



[{^ T -K T d 3 )\ =0. 



Therefore, the Baker-Hausdorff theorem [52] simplifies Eq. (4a) to 

!P(r',i) = e n{r)t e ^ PKat <^ e (r) = e n{r)t e -(^ F ») Kat & e (r) = e n(r)t $ e {r) 



where the last equality again holds because the sum of all internal forces vanishes. 
Therefore, 
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&(r',t) = &(r,t) 

holds, proving that the time-dependent and consequently the stationary distribution 
function & s (T) = lim^oo &(r, t) are translationally invariant. This applies, at least, 
in cases without spontaneous symmetry breaking. Formally, the role of such symme- 
try breaking is to discard some parts of the steady state distribution function and keep 
others (with the choice dependent on initial conditions). The distributions developed 
here discard nothing, and would therefore average over the disjoint symmetry-related 
states of a symmetry-broken system. 

Appreciable simplifications follow from translational invariance for steady-state 
quantities of wavevector-dependent fluctuations: 

/ q (r,t) = e ^ 4 x I(n e iq - r * , 

i 

where e.g. X? = 1 describes density fluctuations g q (i), while X? = 8 a p + 
(l/2)X^-( r ? — rJ)du(\Ti — Tj\)/drf gives the stress tensor element cx a/ 3(q) for 
interactions described by the pair-potential u. Translational invariance in an infinite 
sheared system dictates that averages are independent of identical shifts of all particle 
positions. As the integral over phase space must agree for either integration variables 
r or r', steady-state averages can be non-vanishing for zero wavevector only: 

^{/q(t)> W =/0(7)«q,O. 

The average density n = N/V and the shear stress (7(7) = (c X y)^ /V are impor- 
tant examples. Wavevector-dependent steady-state structure functions under shear 
become anisotropic but remain translationally invariant, so that introduction of a sin- 
gle wavevector suffices. The structure factor built with density fluctuations shall be 
abbreviated by 

It needs to be kept apart from the equilibrium structure factor, denoted by 

Sq = (<5£ q Sgq) , (6b) 

which is obtained by averaging over the particle positions using the equilibrium dis- 
tribution function <P e . It will be one of the hallmarks of a shear molten, yielding glass 
state, that even in the limit of vanishing shear rate both structure factors do not agree: 
Sq(j — > 0) 7^ S q in a shear molten glassy state. 

Translational invariance of sheared systems takes a special form for two-time 
correlation functions, because a shift of the point in coordinate space from r to J" 
gives 

(Sf* e nU 6g k }^ - e -i(k-«t+k-q).a {§f * e ^*^ k )(7) ; 
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t=0 




q(0) 



q(t) - 



Fig. 1. Shear advection of a fluctuation with initial wavevector in 2>direction, q(i = 0) = 
q (1, 0, 0) T , and adverted wavevector at later time q(i>0) = q (1, —7*, 0) T ; from [44]. While 
X x is the wavelength in ^-direction at t = 0, at later time t, the corresponding wavelength A y 
in (negative) -^-direction obeys: \ x /\ y = Ax/ Ay — jt. At all times, q(t) is perpendicular 
to the planes of constant fluctuation amplitude. Note that the magnitude q(t) = qyj 1 + { A /t) 2 
increases with time. Brownian motion, neglected in this sketch, would smear out the fluctua- 
tion. 

while obviously both averages need to agree. Therefore, a fluctuation with wavevec- 
tor q is correlated with a fluctuation of k = q(i) with the advected wavevector 



at the later time t; only then the exponential in the last equation becomes unity; 
fluctuations with other wavevector combinations are decorrelated. The advected 
wavevector's y-component decreases with time as q y (t) = q y — A ft q x , corresponding 
to an (asymptotically) decreasing wavelength, which the shear-advected fluctuation 
exhibits along the y-direction; see Fig. 1. Taking into account this time-dependence 
of the wavelength of fluctuations, a stationary time-dependent correlation function 
characterized by a single wavevector can be denned: 



Application of Eqs. (4,5) is potentially obstructed by the existence of conser- 
vation laws, which may cause a zero eigenvalue of the (adjoint) SO, . The time 
integration in Eqs. (4b, 5) would then not converge at long times. This possible obsta- 
cle when performing memory function integrals, and how to overcome it, is familiar 
from equilibrium Green-Kubo relations [53]. For Brownian particles, only the den- 
sity q is conserved. Yet, density fluctuations do not couple in linear order to the 
shear-induced change of the distribution function [43]. The (equilibrium) average 



q(t) = q — q • nt 



(7) 



(8) 




ftj) =0 



vanishes for all q; at finite q because of translational invariance, and at q = because 
of inversion symmetry. Thus, the projector Q can be introduced 
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Q = 



1 - P , with P = 5Qq ) 



1 



(9) 



q 



It projects any variable into the space perpendicular to linear density fluctuations. 
Introducing it into Eq. (5c) is straight forward, because couplings to linear density 
can not arise in it anyway. One obtains 



" Jo 

The projection step is exact, and also formally redundant at this stage; but it will 
prove useful later on, when approximations are performed. 

2.2.3 

Coupling to structural relaxation 

The generalized Green-Kubo relations, leave us with the problem of how to approxi- 
mate time-dependent correlation functions in Eq. (5). Their physical meaning is that 
at time zero, an equilibrium stress fluctuation arises; the system then evolves under 
internal and shear-driven motion until time t, when its correlation with a fluctuation 
Afq=o is determined. Integrating up these contributions for all times since the start 
of shearing gives the difference of the shear-dependent quantities to the equilibrium 
ones. During the considered time evolution, the projector Q prevents linear couplings 
to the conserved particle density. 

The time dependence and magnitudes of the correlations in Eq. (5) shall now be 
approximated by using the overlaps of both the stress and Afq =0 fluctuations with 
appropriately chosen 'relevant slow fluctuations' . For the dense colloidal dispersions 
of interest, the relevant structural rearrangements are assumed to be density fluctua- 
tions. Because of the projector Q in Eq. (5d), the lowest nonzero order in fluctuation 
amplitudes, which we presume dominant, must then involve pair-products of density 
fluctuations, g p . 

The mode coupling approximation may be summarized as a rule that applies to 
all fluctuation products that exhibit slow structural relaxations but whose variables 
cannot couple linearly to the density. Their time-dependence is approximated as: 



The fluctuating variables are thereby projected onto pair-density fluctuations, whose 
time-dependence follows from that of the transient density correlators <£ q ( t )(t), de- 
fined in Eq. (12) below. These describe the relaxation (caused by shear, interactions 
and Brownian motion) of density fluctuations with equilibrium amplitudes. Higher 
order density averages are factorized into products of these correlators, and the re- 
duced dynamics containing the projector Q is replaced by the full dynamics. The 
entire procedure is written in terms of equilibrium averages, which can then be used 
to compute nonequilibrium steady states via the ITT procedure. The normalization 




(5d) 
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in Eq. (10a) is given by the equilibrium structure factors such that the pair density 
correlator with reduced dynamics, which does not couple linearly to density fluctua- 
tions, becomes approximated to: 

(QtQpQ e Qn ' Qt Q Qk(t) Qp(t) > » N 2 S k S p ^ k (i) * p (t) . (10b) 

This equation can be considered as central approximation of the MCT [37] and 
MCT-ITT approach. While the projection onto density pairs, which is also con- 
tained/implied in Eq. (10a) may be improved upon systematically by including 
higher order density or other fluctuations, see [54, 55] for examples, no systematic 
way to improve upon the breaking of averages in Eq. (10b) has been discovered up 
to now, to the knowledge of the author. 

The mode coupling approximations introduced above can now be applied to the 
exact generalized Green-Kubo relations Eq. (5d). Steady state expectation values are 
approximated by projection onto pair density modes, giving 

m « (f )/v + ± jy ]r ^Mz||p V f g2 ( _ t)(t) (lla) 

with t the time since switch-on of shear. To derive this, the property = ^_k = 
was used; also the restriction k > p when summing over wavevectors was dropped, 
and a factor \ introduced, in order to have unrestricted sums over k. Within Eq. (11a) 
we have already substituted the following explicit result for the equal-time correlator 
of the shear stress with density products: 



.k x k y {-t) _N 

"fcRT k{ ~ t] dk( -* ) '- p( - t) " T 

(lib) 



(vxy Q Qk(-t) Qp(-t)) - N X U V _^ S' k( _ t) 5 k( _ t) -p(-t) - -r d t S q( _ t) 5 k _ p 



It's an exact equality using the equilibrium distribution function and Eq.(6) 

(<r xv Q p k A) = (<r xy Qk gt) = / drv e {- J^F l x y t ) g k g* k 



= J dTf^dfyi) g* gt = ik^iviie^-^ - e"*^-^))) 

i ij 

Equation (11a), as derived via the mode-coupling rule detailed above, contains 
a 'vertex function' V k , describing the coupling of the desired variable / to density 
pairs. This denotes the following quantity, computed using familiar thermodynamic 
equalities 

V( = (g^g* Q Af )/N = (gteu Af )/N-S ( S k + n dSk ^ 1 9{fo) 



dn J V dn 



T 



(lie) 
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In ITT, the slow stress fluctuations in g(t, 7) are approximated by following the 
slow structural rearrangements, encoded in the transient density correlators. The gen- 
eralized modulus becomes, using the approximation Eq. (10a) and the vertex Eq. 
(lib): 

nl+ A,\ keT f d3k ^hkvi-t) S k S H-t) ^2 M nlH x 
ff(*,7) = ~ / 7^)3 -J^Zt] sT~ ^ {t) ' (Ud) 

Summation over wavevectors has been turned into integration in Eq. (lid) consider- 
ing an infinite system. 

The familiar shear modulus of linear response theory describes thermodynamic 
stress fluctuations in equilibrium, and is obtained from Eqs. (5b, lid) by setting 
7 = [1,3,56]. While Eq. (5b) then gives the exact Green-Kubo relation, the approx- 
imation Eq. (lid) turns into the well-studied MCT formula (see Eq. (17) below). For 
finite shear rates, Eq. (lid) describes how affine particle motion causes stress fluctu- 
ations to explore shorter and shorter length scales. There the effective forces, as mea- 
sured by the gradient of the direct correlation function, S' k /S% = nc' k = ndck/dk, 
become smaller, and vanish asympotically, — ► 0; the direct correlation func- 
tion Cfe is connected to the structure factor via the Ornstein-Zernicke equation 
Sk = 1/(1 — n Cfe), where n = N/V is the particle density. Note, that the equilib- 
rium structure suffices to quantify the effective interactions, while shear just pushes 
the fluctuations around on the 'equilibrium energy landscape'. 

While, in the linear response regime, modulus and density correlator are measur- 
able quantities, outside the linear regime, both quantities serve as tools in the ITT 
approach only. The transient correlator and shear modulus provide a route to the sta- 
tionary averages, because they describe the decay of equilibrium fluctuations under 
external shear, and their time integral provides an approximation for the stationary 
distribution function. Determination of the frequency dependent moduli under large 
amplitude oscillatory shear has become possible recently only [57], and requires an 
extension of the present approach to time dependent shear rates in Eq. (3) [58]. 



2.2.4 

Transient density correlator 

In ITT, the evolution towards the stationary distribution at infinite times is approx- 
imated by following the slow structural rearrangements, encoded in the transient 
density correlator ^ q (£). It is defined by [42,43] 

*q(*)= j^r<^e fltt «y eq(t) )W=°). (12) 

It describes the fate of an equilibrium density fluctuation with wavevector q, where 
g q = J2f=i e lq rj , under the combined effect of internal forces, Brownian motion 
and shearing. Note that because of the appearance of fy e in Eq. (4), the average in 
Eq. (12) can be evaluated with the equilibrium canonical distribution function, while 
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the dynamical evolution contains Brownian motion and shear advection. The nor- 
malization is given by S q the equilibrium structure factor [3, 45] for wavevector 
modulus q = |q|. The advected wavevector from Eq. (7) enters in Eq. (12). The 
time-dependence in q(i) results from the affine particle motion with the shear flow 
of the solvent. Again, irrespective of the use of ^ e in Eq. (12), or \P S in Eq. (8), in 
both cases translational invariance under shear dictates that at a time t later, the den- 
sity fluctuation <5g* has a nonvanishing overlap only with the advected fluctuation 
(5g q ( t ) . Figure 1 again applies, where a non-decorrelating fluctuation is sketched un- 
der shear. In the case of vanishing Brownian motion, viz. Do = in Eq. (3b), we 
find ^q(t) = 1, because the advected wavevector takes account of simple affine par- 
ticle motion. The relaxation of <£ q (t) thus heralds decay of structural correlations by 
Brownian motion, affected by shear. 

2.2.5 

Zwanzig-Mori equations of motion 

Structural rearrangements of the dispersion affected by Brownian motion is encoded 
in the transient density correlator. Shear induced affine motion, viz. the case Dq = 0, 
is not sufficient to cause ^k(i) to decay. Brownian motion of the quiescent correlator 
$k~°\t) leads at high densities to a slow structural process which arrests at long 
times in (metastable) glass states. Thus the combination of structural relaxation and 
shear is interesting. The interplay between intrinsic structural motion and shearing 
in $k(£) is captured by (i) first a formally exact Zwanzig-Mori type equation of 
motion, and (ii) second a mode coupling factorisation in the memory function built 
with longitudinal stress fluctuations [A2-AA]. The equation of motion for the transient 
density correlators is 

9t* q (t) + r q (t) j«P q (t) + J dt' m q (t,t') b\> ^q(i')} - , (13) 

where the initial decay rate 2"q(t) = Do q 2 (t)/S q ( t ) generalizes the familiar result 
from linear response theory to advected wave vectors; it contains Taylor dispersion 
mentioned in the introduction, and describes the short time behavior, <Pq(t — > 0) — > 

i - r q (o) t + . . . . 

2.2.6 

Mode-coupling closure 

The memory equation contains fluctuating stresses and similarly like g(t, 7) in Eq. 
(lid), is calculated in mode coupling approximation using Eq. (10a) giving: 

"»,(*, *0 = ^ 2 W*' *') *k(f )(* - t') # p(t<) (t - t') , (14a) 

k 

where we abbreviated p = q k. The vertex generalizes the expression in the qui- 
escent case, see Eq. (18c) below, and depends on two times capturing that shearing 
decorrelates stress fluctuations [42-44] 
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wmo = Jq( ;^ ( y o v qkp( t) v qkP (f) , 

V q kp(t) = q(t) • ( k(i) nc fc( t) + p(t) nc p(t) ) . (14b) 

With shear, wavevectors in Eq. (14b) are advected according to Eq. (7). 

The summarized MCT-ITT equations form a closed set of equations determining 
rheological properties of a sheared dispersion from equilibrium structural input [42- 
44]. Only the static structure factor S q is required to predict (i) the time dependent 
shear modulus within linear response, g lr (t) — g(t, 7 = 0), and (ii) the stationary 
stress (7(7) from Eq. (5a). 



2.3 

A microscopic model: Brownian hard spheres 

In the microscopic ITT approach, the rheology is determined from the equilibrium 
structure factor S q alone. This holds at low enough frequencies and shear rates, and 
excludes a single time scale, to be denoted by the parameter to in Eq. (22b), which 
needs to be found by matching to the short time dynamics. This prediction has as 
consequence that the moduli and flow curves should be a function only of the ther- 
modynamic parameters characterizing the present system, viz. its structure factor. 
Because the structure factor for simple fluids far away from demixing and other 
phase separation regions can be mapped onto the one of hard spheres, the system 
of hard spheres plays a special role in the MCT-ITT approach. It provides the most 
simple microscopic model where slow strucutral dynamics can be studied. Moreover, 
other experimental systems can be mapped onto it by chosing an effective packing 
fraction = {A-K/2i)nF? H and particle radius so that the structure factors agree. 

The claim that the rheology follows from S q is supported if the rheological prop- 
erties of a dispersion only depend on the effective packing fraction, if particle size 
is taken account of properly. Obviously, appropriate scales for frequency, shear rate 
and stress magnitudes need to be chosen to observe this; see Sect. 6.2. The depen- 
dence of the rheology (via the vertices) on S q suggests that ksT sets the energy scale 
as long as repulsive interactions dominate the local packing. The length scale is set 
by the average particle separation, which can be taken to scale with Rh- The time 
scale of the glassy rheology within ITT is given by to, which should scale with the 
measured dilute diffusion coefficient Do- Thus the rescaling of the rheological data 
can be done with measured parameters alone. 

Because the hard sphere system thus provides the most simple system to test and 
explore MCT-ITT, numerical calculations only for this model will be reviewed in the 
present overview. Input for the structure factor is required, which, for simplicity, will 
be taken from the analytical Percus-Yervick (PY) approximation [2, 3]. Straightfor- 
ward discretization of the wavevector integrals will be performed as discussed below, 
and in detail in the quoted original papers. 
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2.4 

Accounting for hydrodynamic interactions 

Solvent-particle interactions (viz. the HI) act instantaneously if the particle mi- 
crostructure differs from the equilibrium one, but do not themselves determine the 
equilibrium structure [3,45]. If one assumes that glassy arrest is connected with the 
ability of the system to explore its configuration space and to approach its equilib- 
rium structure, then it appears natural to assume that the solvent particle interactions 
are characterized by a finite time scale thi. And that they do not shift the glass tran- 
sition nor affect the frozen glassy structure. HI would thus only lead to an increase 
of the high frequency viscosity above the solvent value; this value shall be denoted 
as r/oc: 

ff(t,7)->s(t,7)+»?oo5(t-0+). < 15a > 

The parameter ijoo would thus characterize a short-time, high frequency viscosity and 
model viscous processes which require no structural relaxation. It can be measured 
from the high frequency dissipation 

G"(lo -» oo) = r/oo w . (15b) 

For identical reasoning, also the short time diffusion in the collective (and single 
particle) motion will be affected by HI. The most simple approximation is to adjust 
the initial decay rate 

r^(t) = D s q 2 {t)/S q{t) , (15c) 

where the collective short time diffusion coefficient D s accounts for HI and other 
(almost) instantaneous effects which affect the short time motion, and which are not 
explicitly included in the MCT-ITT approach. 

The naive picture sketched here, is not correct for a number of reasons. It is well 
known that for hard spheres without HI the quiescent shear modulus diverges for 
short times, g lr > HSnoHI (£ — > 0) <~ t -1 / 2 . Lubrication forces, which keep the particles 
apart, eliminate this divergence and render g h ' m (t — > 0) finite [59]. Thus, the simple 
separation of the modulus into HI and potential part is not possible for short times, 
at least for particles with a hard core. Moreover, comparison of simulations without 
and with HI has shown that the increase of (770 — ??oo)/f7oo depends somewhat on HI, 
and thus not just on the potential interactions as implied. 

Nevertheless the sketched picture provides the most basic view of a glass tran- 
sition in colloidal suspensions, connecting it with the increase of the structural re- 
laxation time r. Increased density or interactions cause a slowing down of particle 
rearrangements which leave the HI relatively unaffected, as these solvent mediated 
forces act on all time scales. Potential forces dominate the slowest particle rearrange- 
ments because vitrification corresponds to the limit where they actually prevent the 
final relaxation of the microstructure. The structural relaxation time r diverges at the 
glass transition, while r HI stays finite. Thus close to arrest a time scale separation is 
possible, t 3> t hi . 
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2.5 

Comparison with other MCT inspired approaches to sheared fluids 

The MCT-ITT approach aims at describing the steady state properties of a concen- 
trated dispersion under shear. Stationary averages are its major output, obtained via 
the integration through transients procedure from (approximate) transient fluctuation 
functions, whose strength is the equilibrium one, and whose dynamics originates 
from the competition between Brownian motion and shear induced decorrelation. In 
this respect, the MCT-ITT approach differs from the interesting recent generalization 
of MCT to sheared systems by Miyazaki, Reichman and coworkers [39,40]. These 
authors considered the stationary but time-dependent fluctuations around the steady 
state, whose amplitude is the stationary correlation function, e.g. in the case of den- 
sity fluctuations, it is the distorted structure factor S q (7) Eq. (6a). In the approach 
by Miyazaki et al. this structure factor is an input-quantity required to calculate the 
dynamics, while it is an output quantity, calculated in MCT-ITT from the equilibrium 
S q of Eq. (6b). Likewise, the stationary stress as function of shear rate, viz. the flow 
curve (7(7), is a quantity calculated in MCT-ITT, albeit using mode coupling approx- 
imations, while in the approach of Refs. [39,40] additional ad-hoc approximations 
beyond the mode coupling approximation are required to access <r(j). Thus, while 
the scenario of an non-equilibrium transition between a shear-thinning fluid and a 
shear molten glass, characterized by universal aspects in e.g. (7(7) — see the discus- 
sion in Sect. 4 — forms the core of the MCT-ITT results, this scenario can not be 
directly addressed based on Refs. [39,40]. 

Because the recent experiments and simulations reviewed here concentrated on 
the universal aspects of the novel non-equilibrium transition, focus will be laid on the 
MCT-ITT approach. Reassuringly, however, many similarities between the MCT-ITT 
equations and the results by Miyazaki and Reichman exist, even though these authors 
used a different, field theoretic approach to derive their results. This supports the ro- 
bustness of the mechanism of shear-advection in Eq. (7) entering the MCT vertices in 
Eqs. (1 ld,14), which were derived independently in Refs. [39,40] and Refs. [42-A4-] 
from quite different theoretical routes. This mechanism had been known from ear- 
lier work on the dynamics of critical fluctuations in sheared systems close to phase 
transition points [60], on current fluctuations in simple liquids [61], and on inco- 
herent density fluctuations in dilute solutions [62]. Different possibilities also exist 
to include shear into MCT-inspired approaches, especially the one worked out by 
Schweizer and coworkers including strain into an effective free energy [41]. This ap- 
proach does not recover the (idealized) MCT results reviewed below but starts from 
the extended MCT where no true glass transition exists and describes a crossover 
scenario without e.g. a true dynamic yield stress as discussed below. 

3 

Microscopic results in linear response regime 

Before turning to the properties of the stationary non-equilibrium states under shear, 
it is useful to investigate the quiescent dispersion close to vitrification. Consensus on 
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the ultimate mechanism causing glassy arrest may yet be absent, yet, the so-called 
'cage effect' has lead to a number of fruitful insights into glass formation in dense 
colloidal dispersions. For example, it was extended to particles with a short ranged 
attraction leading to at first surprising predictions [63-66]. 




Fig. 2. Positions of hard disks in two dimensions from a Monte Carlo simulation at a density 
close to freezing; courtesy of Th. Franosch. A particle and its shell of neighbours is highlighted 
by different colors/shadings. 

Figure 2 shows a section of the cell of a Monte Carlo simulation of hard disks 
moving in two dimensions (for simplicity of visualization). The density is just below 
freezing and the sample was carefully equilibrated. Only 100 disks were simulated, 
so that finite size effects cannot be ruled out. Picking out a disk, it is surrounded 
by a shell of on average 6 neighbours (in two dimension, of 12 neighbours in three 
dimension), which hinder its free motion. In order for the central particle to dif- 
fuse at long times, it needs to escape the shell of neighbours. In order for a gap in 
this shell to open at higher concentrations, the neighbours have to be able to move 
somewhat themselves. Yet, each neighour is hindered by its own shell of neighbours, 
to which the originally picked particle belongs. Thus one can expect a cooperative 
feedback mechanism that with increasing density or particle interactions particle re- 
arrangements take more and more time. It appears natural, that consequently stress 
fluctuations also slow down and the system becomes viscoelastic. 

MCT appears to capture the cage-effect in supercooled liquids and predicts that 
it dominates the slow relaxation of structural correlations close to the glass transition 
[2, 37]. Density fluctuations play an important role because they are well suited to 
describe the structure of the particle system and its relaxation. Moreover, stresses that 
decay slowly because of the slow particle rearrangements, MCT argues, also can be 
approximated by density fluctuations using effective potentials. Density fluctuations 
not at large wavelengths, but for wavelengths corresponding to the average particle 
distance turn out to be the dominant ones. In agreement with the picture of the caging 
of particles by structural correlations, the MCT glass transition is independent on 
whether the particles move ballistically in between interactions with their neighbors 
(say collisions for hard spheres) or by diffusion. Structural arrest happens whenever 
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the static density correlations for wavelengths around the average particle distance 
are strong enough. The arrest of structural correlations entails an increase in the 
viscosity of the dispersion connected to the existence of a slow Maxwell-process 
in the shear moduli. While the MCT solutions for density fluctuations have been 
thoroughly reviewed, the viscoelastic spectra have not been presented in such detail. 



3.1 

Shear moduli close to the glass transition 
3.1.1 

MCT equations and results for hard spheres 

The loss and storage moduli of small amplitude oscillatory shear measurements [1,3] 
follow from Eq. (5b) in the linear response case at 7 = 0: 

/•OC 

G'(w)+iG"(u)=iu / dte- lwt g lr (t) . (16a) 
Jo 

Here, the shear modulus in the linear response regime is, again like the transient one 
in Eq. (5b), obtained from equilibrium averaging: 

- 7^7 ( v*y ^ v xy )<-i=V , (16b) 

yet, differently from the transient one, the equilibrium one contains the equilibrium 
SO J? e , which characterizes the quiescent system: 

The linear response modulus thus quantifies the small stress fluctuations, which are 
excited by thermally, and relax because of Brownian motion. 

Predictions of the (idealized) MCT equations for the potential part of the equi- 
librium, time-dependent shear modulus g k (t) of hard spheres for various packing 
fractions <f> are shown in Fig. 3 and calculated from the limit of Eq. (lid) for vanish- 
ing shear rate: 

The normalized density fluctuation functions are calculated self-consistently within 
MCT from the Eqs. (13,14) at vanishing shear rate [2, 37], which turn into the quies- 
cent MCT equations: 

d t $ q {t) + r q |<P,(t) + J dt' m q {t - t') d v = , (18a) 
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Fig. 3. Equilibrium potential or linear response shear modulus G(t) = g h {t) (in units of 
ksT/Rfj) for Brownian hard spheres with radius a = Rh for packing fractions close to 
vitrification at <f> c ; from Ref. [67]. Densities are measured by the separation parameter e = 
(<f> — (f>c)/<t>c = ilO - ' 71 '^ 3 , and labels denote the value n. Positive values belong to glass 
(e > 0), negative to fluid states (e < 0); the label c gives the transition. The inset shows 
a subset of the curves on a linear time axis; the increase of g (t) for short times cannot be 
resolved. 



where the initial decay rate r q = D s q 2 /S q describes diffusion with a short-time 
diffusion coefficient D s differing from the Dq because of HI; D s = Do will be taken 
for exemplary calculations, while D s ^ Dq is required for analyzing experimental 
data. The memory kernel becomes (again with abbreviation p = q k) 

m * (*) = *W ^ (*) ^ W ' < 1 8b > 

k 
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qkp 




(q • [k nc k + p ncp}) 



(18c) 



Packing fractions are conveniently measured in relative separations e = ((f)— <j> c ) / <f> c 



Fig. 4. Storage part of the shear modulus G'(lo) corresponding to Fig. 3. The inset shows 
storage and loss moduli (only for fluid states) for a number of densities. 

to the glass transition point, which for this model of hard spheres lies at cj) c = 0.516 
[37,68]. Note that this result depends on the static structure factor S(q) , which 
is taken from Percus-Yevick theory, and that the experimentally determined value 
^expt. _ o 5g jj es somewhat higher [13, 14]. The wavevector integrals were dis- 
cretized using M = 100 wavevectors chosen from fc m i„ = 0.1/ Rh up to k niax — 
19.9/Rh with separation Ak = 0.2/R H for Figs. 3 to 5, or using M = 600 
wavevectors chosen from fc m j n = 0.05/Rh up to k max = 59.95/Rh with sepa- 
ration Ak = 0.1/ Rh for Figs. 6 and 7, and in Sect. 6.2 in Fig. 23. Time was dis- 
cretized with initial step-width dt = 2 10~ 7 R H / D s , which was doubled each time 
after 400 steps. Slightly different discretizations in time and wavevector of the MCT 
equations were used in Sects. 3.2. 1 and 4, causing only small quantitative differences 
whose discussion goes beyond the present review. The quiescent density correlators 
<£ g (i) corresponding to the following linear response moduli have thoroughly been 
discussed in Ref. [68]. 

For low packing fractions, or large negative separations, the modulus decays 
quickly on a time-scale set by the short-time diffusion of well separated particles. 
The strength of the modulus increases strongly at these low densities, and its behav- 
ior at short times presumably depends sensitively on the details of hydrodynamic 
and potential interactions; thus Fig. 3 is not continued to small times, where the em- 
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Fig. 5. Loss part of the shear modulus G"(ui) corresponding to Fig. 3. The inset shows the 
same data in a double logarithmic representation. 

ployed model (taken from Ref. [68,69]) is too crude 3 . Approaching the glass tran- 
sition from below, e /" 0, little changes in g lr (t) at short times, because the absolute 
change in density becomes small. Yet, at long times a process in g lr (t) becomes pro- 
gressively slower upon taking e to zero. It can be considered the MCT analog of 
the phenomenological Maxwell-process. MCT finds that it depends on the equilib- 
rium structural correlations only, while HI and other short time effects only shift its 
overall time scale. Importantly, this overall time scale applies to the slow process in 
coherent and incoherent density fluctuations as well as in the stress fluctuations [70]. 
This holds even though e.g. HI are known to affect short time diffusion coefficients 
and high frequency viscosities differently. Upon crossing the glass transition, a part 
of the relaxation freezes out and the amplitude Goo of the Maxwell-process does 
not decay; the modulus for long times approaches the elastic constant of the glass 
g lr (t — ► oo) — > Goo > 0. Entering deeper into the glassy phase the elastic constant 
increase quickly with packing fraction. 

The corresponding storage G'(ui) and loss G"(ui) moduli are shown as functions 
of frequency in Figs. 4 and 5, respectively. The slow Maxwell-process appears as a 
shoulder in G' which extends down to lower and lower frequencies when approach- 
ing glassy arrest, and reaches to zero frequency in the glass, G'(u> = 0) = Goo. 
The slow process shows up as a peak in G" which in parallel motion (see inset of 
Fig. 4) shifts to lower frequencies when e /" 0. Including hydrodynamic interactions 

A The MCT shear modulus at short times depends sensitively on the large cut-off fc max for 
hard spheres [56], g(t,-y = 0) = (n 2 k B T/60n 2 ) /'"" dfcfc 4 (4) 2 Sl$\{t) gives the 

qualitatively correct [59,71] short time g k (t — > 0) ~ t -1 ^ 2 , or high frequency divergence 

G'(ui » Dq/Rh) ~ y/uj only for fc max — ► oo. 
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into the calculation by adjusting 77^ would affect the frequency dependent moduli 
at higher frequencies only. For the range of smaller frequencies which is of interest 
here, only a small correction would arise. 



3.1.2 

Comparison with experiments 

Recently, it has been demonstrated that suspensions of thermosensitive particles 
present excellent model systems for measuring the viscoelasticity of dense concen- 
trations. The particles consist of a solid core of polystyrene onto which a thermosen- 
sitive network of poly(N-isopropylacrylamide) (PNIPAM) is attached [30, 31]. The 
PNiPAM shell of these particles swells when immersed in cold water (10 - 15°C). 
Water gets expelled at higher temperatures leading to a considerable shrinking. Thus, 
for a given number density the effective volume fraction e ff can be adjusted within 
wide limits by adjusting the temperature. Senff et cd. (1999) were the first to demon- 
strate the use of these particles as model system for studying the dynamics in concen- 
trated suspensions [23,24]. The advantage of these systems over the classical hard 
sphere systems are that dense suspensions can be generated in situ without shear and 
mechanical deformation. The previous history of the sample can be erased by raising 
the temperature and thus lowering the volume fraction to the fluid regime. 




Fig. 6. The reduced storage (diamonds and solid lines) and loss (squares and broken lines) 
modulus for a fluid state at effective packing fraction cf> e f[ = 0.540; from Ref. [31]. The 
vertical bars mark the minimal rescaled frequency above which the influence of crystallisation 
can be neglected. Parameters in the MCT calculation given as blue lines: e = —0.01, — 
0.15, and 7]^ = 0.3 ksT / (DqRh); moduli scale factor c y = 1.4. For the other lines see 
Ref. [31]. 
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Fig. 7. The moduli for a glass state at effective packing fraction tf> e g = 0.622; storage (di- 
amonds and solid lines) and loss (squares and broken lines) moduli from [31]. MCT fits are 
shown as blue lines with parameters: e = 0.03, = 0.08, and r^oo = 0.3 fcsT /(DoRh); 
moduli scale factor c y — 1.4 . For the other lines see Ref. [31]. 



Frequency dependent moduli were measured spanning a wide density and fre- 
quency range by combining different techniques. The moduli exhibit a qualitative 
change when increasing the effective packing fraction from around 50% to above 
60%. For lower densities (see Fig. 6), the spectra G"(lo) exhibit a broad peak or 
shoulder, which corresponds to the final or a-relaxation. Its peak position (or alter- 
natively the crossing of the moduli, G = G") is roughly given by lot = 1 . These 
properties characterize a viscoelastic fluid. For higher density, see Fig. 7, the storage 
modulus exhibits an elastic plateau at low frequencies. The loss modulus drops far 
below the elastic one. These observations characterize a soft solid 4 . 

The linear response moduli are affected by the presence of small crystallites. At 
low frequencies, G'(lo) and G"{lo) increase above the behavior expected for a solu- 
tion (G'(lj — ► 0) — * r/o uj and G"(oj — ► 0) — > clj 2 ) even at low density, and exhibit 
elastic contributions (apparent from G'(ui) > G"(ui)); see Fig. 6. This effect tracks 
the crystallisation of the system during the measurement after a strong preshearing at 
7 = 100 s _1 . Only data can be considered which were collected before the crystalli- 
sation time; they lie to the right of the vertical bar in Fig. 6. While this experimental 
restriction limits more detailed studies of the shapes of the spectra close to the glass 
transition, the use of a system with a rather narrow size distribution provides the 
quantitatively closest comparison with MCT calculations for a monodisperse hard 

4 The loss modulus rises again at very low frequencies, which may indicate that the colloidal 
glass at this density is metastable and may have a finite lifetime (an ultra-slow process is 
discussed in Ref. [31]). 
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sphere system. Especially the magnitude of the stresses and the effective densities 
can be investigated quantitatively. 

Included in figures 6 and 7 are calculations using the microscopic MCT given by 
Eqs. (16) to (18) evaluated for hard spheres in PY approximation. The only a priori 
unknown, adjustable parameter is a frequency or time scale, which was adjusted by 
varying the short time diffusion coefficient D s appearing in the initial decay rate in 
Eq. (18a). Values for D s /Dq are reported in the captions. The viscous contribution 
to the stress is mimicked by including r/oo like in Eq. (15); it can directly be mea- 
sured at the highest frequencies. Gratifyingly, the stress values computed from the 
microscopic approach are close to the measured ones; they are too small by 40% 
only, which may arise from the approximate structure factors entering the MCT cal- 
culation; the Percus-Yevick approximation was used here [3]. In order to compare 
the shapes of the moduli the MCT calculations were scaled up by a factor c y = 1.4 
in Figs. 6 and 7. Microscopic MCT also does not hit the correct value for the glass 
transition point [2, 37]. It finds <^ CT = 0.516, while experiments give 0° xp w 0.58. 
Thus, when comparing, the relative separation from the respective transition point 
needs to be adjusted as, obviously, the spectra depend sensitively on the distance to 
the glass transition; the fitted values of the separation parameter e are included in the 
captions. 

Overall, the semi-quantitative agreement between the linear viscoelastic spectra 
and first-principles MCT calculations is very promising. Yet, crystallization effects 
in the data prevent a closer look, which will be given in Section 6.2, where data from 
a more polydisperse sample are discussed. 

3.2 

Distorted structure factor 
3.2.1 

Linear order in 7 

The stationary structural correlations of a dense fluid of spherical particles undergo- 
ing Brownian motion, neglecting hydrodynamic interactions, change with shear rate 
7 in response to a steady shear flow. In linear order, the structure is distorted only in 
the plane of the flow, while already in second order in 7, the structure factor changes 
under shear also for wavevectors lying in the plane perpendicular to the flow. Consis- 
tent with previous theories, MCT-ITT finds regular expansion coefficients in linear 
and quadratic order in 7 for fluid (ergodic) suspensions [72]. For the steady state 
structure factor 6^(7) of density fluctuations under shearing in plain Couette flow 
defined in Eq. (6a), the change from the equilibrium one in linear order in shear rate 
is given by the following ITT-approximation: 

5 9 ( 7 ) = S q +7 {^fs> q Jit *2(t)} + G(Pe 2 ) . (19) 

This relation follows from Eq. (11a) in the limit of small 7, where the quiescent 
density correlators can be taken from quiescent MCT in Eq. (18). 



Nonlinear rheology 



29 



3.2.2 

Comparison with simulations 

Figure 8 shows the contribution SSq to the distorted structure factor in leading linear 
order in 7 for packing fractions <j> = 0.36, 0.44 and 0.46. Data taken from Brown- 
ian dynamics simulations at (p = 0.43 and 0.5 from Ref. [73] are also included. In 
both cases the data was divided by a factor jq x q y /q 2 which is the origin of the triv- 
ial anisotropy in the leading linear order. The distortion SSq^ of the microstructure 
grows strongly with <f>, because of the approach to the glass transition. The SSq* 1 is 
proportional to the a-relaxation time r, as proven in the left inset of Fig. 8. Here, r 
is estimated from ^ 9p (t = r) = 0.1, where q p denotes the position of the primary 
peak in S q . Rescaling the data with Pe, collapses the curves at different distances to 
the glass transition. The strongest shear-dependence occurs for the direction of the 
extensional component of the flow, q x = q y . Here, the mesoscale order of the disper- 
sion grows; the peak in SSq^ increases and sharpens. The ITT results qualitatively 
agree with the simulations in these aspects [73]. 

The most important finding of Fig. 8 concerns the magnitude of the distortion of 
the microstructure, and the dimensionless parameter measuring the effect of shear 
relative to the intrinsic particle motion. This topic can already be discussed using 
the linear order result, and is not affected by considerations of hydrodynamic in- 
teractions, as can be glanced from comparing Brownian dynamics simulations [73] 
and experiments on dissolved particles [74]. In previous theories, shear rate effects 
enter when the bare Peclet number Peo becomes non-negligible. In the present ITT 
approach the dressed Peclet/ Weissenberg number Pe= 7T governs shear effects; 
here, r is the (final) structural relaxation time. Shear flow competes with structural 
rearrangements that become arbitrarily slow compared to diffusion of dilute particles 
when approaching the glass transition. The distorted microstructure results from the 
competition between shear flow and cooperative structural rearrangements. It is thus 
no surprise that previous theories using Pe , which is characteristic for dilute fluids 
or strong flows, had severely underestimated the magnitude of shear distortions in 
hard sphere suspensions for higher packing fractions; Refs. [73, 75] report an un- 
derestimate by roughly a decade at (j> = 0.50. The ITT approach actually predicts 
a divergence of lim^^ (5'q^' — S q )/j for density approaching the glass transition 
at (j> c :=y 0.58. And for (idealized) glass states, where t = 00 holds in MCT follow- 
ing Maxwell's phenomenology, the stationary structure factor becomes non-analytic, 
and differs from the equilibrium one even for 7^0. The distortion 8Sq^ thus 
qualitatively behaves like the stress, which goes to zero linear in 7 in the fluid, but 
approaches a yield stress <r + in the glass for 7^0. 

The reassuring agreement of ITT results on Sq ^ with the data from simulations 
and experiments shows that in the ITT approach the correct expansion parameter Pe 
has been identified. This can be taken as support for the ITT-strategy to connect the 
non-linear rheology of dense dispersions with the structural relaxation studied at the 
glass transition. 



30 



Matthias Fuchs 




Fig. 8. Contribution SS^ linear in shear rate to the distorted structure factor Sq(j) normal- 
ized to 'yq x q y /q 2 (black solid lines); from Ref. [72], here D denotes the hard sphere diameter. 
Decreasing the relative separations from the critical point as e = —0.2, —0.15, —0.1, the 
magnitude of 6Sq increases. The dashed-dotted red curves are Brownian dynamics simu- 
lation data from Ref. [73] using the same normalization at e = —0.259, —0.138. The right 
inset shows the unnormalized Sq along the extensional axis q x = q y = g/v^ate = -0.1 
and, from bottom to top, Pe = yr = (green-solid), 1/8, 1/4 and 1/2 (all black-dashed), 
where <& qp (t = t) = 0.1 defines r. Pe/Peo = 1.66 holds at this e. The left inset shows the 
data of the main figure rescaled with the dressed Peclet number, 5S^ ,wiso ' /{Pe q x q y /q 2 ); 
for e — —0.1, —0.05, and —0.01 (with increasing peak height) the values Pe/Peo = 1.66, 
Pe/Peo = 8.06, and Pe/Pe = 419 are used. 



4 

Universal aspects of the glass transition in steady shear 

The summarized microscopic MCT-ITT equations contain a non-equilibrium transi- 
tion between a shear thinning fluid and a shear-molten glassy state; it is the central 
novel transition found in MCT-ITT [42]. Close to the transition, (rather) universal 
predictions can be made about the non-linear dispersion rheology and the steady 
state properties. Following Refs. [31,76], the central predictions are introduced in 
this section and summarized in the overview figure 9; the following results sections 
contain more examples. Figure 9 is obtained from the schematic model which is also 
often used to analyse data, and which is introduced in the following section 5.2 . 

A dimensionless separation parameter e measures the distance to the transition 
which lies at e = 0. A fluid state (e < 0) possesses a (Newtonian) viscosity, 770 (s < 
0) = lim^o cr(7)/7, and shows shear-thinning upon increasing 7. Via the relation 
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Fig. 9. Overview of the properties of the F\l -model characteristic for the transition between 
fluid and yielding glass; from Ref. [31]. The upper panel shows numerically obtained transient 
correlators $ (t) for e = 0.01 (black curves), e = (red), e = —0.005 (green), and e = 
-0.01 (blue). The shear rates are \ = (thick solid lines), \ = 10" 6 (dotted 

lines), and |7/-T| = 10 -2 (dashed lines). For the glass state at e = 0.01 (black), \ = 
10 -8 (dashed-dotted-line) is also included. All curves were calculated with 7 C = 0.1 and 
f]oo = 0. The thin solid lines give the factorization result Eq. (20) with scaling functions Q for 
\i/r\ — 10~ 6 ; label a marks the critical law (22b), and label b marks the von Schweidler- 
law (23). The critical glass form factor f c is indicated. The inset shows the flow curves for 
the same values for e. The thin black bar shows the yield stress for e = 0. The lower 
panel shows the viscoelastic storage (solid line) and loss (broken line) modulus for the same 
values of e. The thin green lines are the Fourier-transformed factorization result Eq. (20) with 
scaling function Q taken from the upper panel for e = —0.005. The dashed-dotted lines show 
the fit formula Eq. (24) for the spectrum in the minimum-region with G m m/i> CT = 0.0262, 
LOmin/T = 0.000457 at e = -0.005 (green) and G min /u CT = 0.0370, u mia /r = 0.00105 
at e — —0.01 (blue). The elastic constant at the transition is marked also, while the high 
frequency asymptote G'oc = G'(u) — > oo) is not labeled explicitly. 
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?7o = lim^^o G"(u)/ui, the Newtonian viscosity can also be taken from the linear 
response loss modulus at low frequencies, where G"(u>) dominates over the storage 
modulus. The latter varies like G'(io — > 0) <~ lo 1 . A glass (e > 0), in the absence of 
flow, possesses an elastic constant G^, which can be measured in the elastic shear 
modulus G'(u)) in the limit of low frequencies, G'(co — > 0,e > 0) — > Goo(e). 
Here the storage modulus dominates over the loss one, which drops like G" (u> — > 
0) <~ lo. The high frequency modulus G'^ = G'(uj — > oo) is characteristic of the 
particle interactions, see Footnote 4, and exists, except for the case of hard sphere 
interactions without HI, in fluid and solid states. The dissipation at high frequencies 
G"(u) — > oo) — ► rjoaLj also shows no anomaly at the glass transition and depends 
strongly on HI and solvent friction. 

Enforcing steady shear flow melts the glass. The stationary stress of the shear- 
molten glass always exceeds a (dynamic) yield stress. For decreasing shear rate, the 
viscosity increases like I/7, and the stress levels off onto the yield-stress plateau, 
£7(7 -> 0,e > 0) -> a+{e). 

Close to the transition, the zero-shear viscosity 770, the elastic constant G M , and 
the yield stress a + show universal anomalies as functions of the distance to the tran- 
sition. 




Fig. 1 0. Glass form factors /, as function of wavevector q in a colloidal glass of hard spheres 
for packing fractions cf> as labeled. Data obtained by van Megen and coworkers by dynamic 
light scattering are qualitatively compared to MCT computations using the PY-S 9 at <f> values 
chosen ad hoc to match the experimental data; from [12]. The PY structure factor at the glass 
transition density <f> = 0.58 is shown as broken line, rescaled by a factor 1/10 . 



The described results follow from the stability analysis of Eqs. (13,14) around 
an arrested, glassy structure f q of the transient correlator [42,76]. Considering the 
time window where *P q (i) is close to arrest at f q , and taking all control parameters 
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like density, temperature, etc. to be close to the values at the transition, the stability 
analysis yields the 'factorization' between spatial and temporal dependences 

*q(t) = f c q + K Q{ t/t ,e, t*o ) + ■ ■ ■ , (20) 

where the (isotropic) glass form factor /° and critical amplitude h q describe the 
spatial properties of the metastable glassy state. The critical glass form factor f° 
gives the long-lived component of density fluctuations right at the transition e = 0, 
while <Pq(t — > oo,e > 0, 7 = 0) = f q > characterizes states even deep in the 
glass with fq obeying [37]: 



(21) 



with V q k P from Eq. (18c) as follows from Eqs. (13,14) asymptotically in the limit of 
vanishing shear rate [76]. Figure 10 shows dynamic light scattering data for the glass 
form factors at a number of densities in PMMA hard sphere colloids, comparing them 
to solutions of Eq. (21) evaluated for hard spheres using the PY structure factor; it 
is included in Fig. 10 for the packing fraction <p c = 0.58 of the experimental glass 
transition. The glass form factor varies with the average particle separation and in 
phase with the structure factor. Both, f q and h q thus describe local correlations, the 
so-called 'cage-effect', and can be taken as constants independent on shear rate and 
density, as they are evaluated from the vertices in Eq. (14) at the transition point 
( £ = 0,7 = 0). 

All time-dependence and (sensitive) dependence on the external control parame- 
ters is contained in the function Q, which often is called '/3-correlator' and obeys the 
non-linear stability equation [42, 44, 76] 

e - C ™ {itf + \g\t) = j t Jdt' Q{t - t') G(t') , (22a) 

with initial condition 

0(t->O) -> (t/t y a . (22b) 

The two parameters A and in Eq. (22a) are determined by the static structure 
factor at the transition point, and take values around A w 0.73 and ~ 0.7 for 
the PY S q for hard spheres. The transition point then lies at packing fraction <j> c — 
^-n c R? H « 0.52 (index c for critical), and the separation parameter measures the 
relative distance, e = C e with e = (<j)—(f> c )/(j) c and C sa 1.3. The 'critical' exponent 
a is given by the exponent parameter A via A = _T(1 — a) 2 /r(l — 2a), as had been 
found in quiescent MCT [2, 37]. 

The time scale to in Eq. (22b) provides the means to match the function Q (t) 
to the microscopic, short-time dynamics. The Eqs. (13,14) contain a simplified 
description of the short time dynamics in colloidal dispersions via the initial de- 
cay rate -T q (f). From this model for the short-time dynamics, the time scale t « 
1.6 1Q~ 2 R 2 H / Do is obtained. Solvent mediated effects on the short time dynamics 
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are well known and are neglected in r q (i) in Eq. (13). The most simple minded 
approximation to account for HI is given in Eq. (15). It only shifts the value of t - 
Within the ITT approach, this finding holds more generally. Even if HI lead to more 
substantial changes of Eq. (13), all of the mentioned universal predictions would re- 
main true, as long as HI do not affect the mode coupling vertex in Eq. (14). Like in 
the quiescent MCT [70], in MCT-ITT hydrodynamic interactions can thus be incor- 
porated into the theory of the glass transition under shear, and amount to a rescaling 
of the matching time t , only. 

The parameters e, A and in Eq. (22a) can be determined from the equilib- 
rium structure factor S q at or close to the transition, and, together with to and the 
shear rate 7 they capture the essence of the rheological anomalies in dense disper- 
sions. A divergent viscosity follows from the prediction of a strongly increasing final 
relaxation time in Q in the quiescent fluid phase [2, 37]: 

Q(t -» 00, e < 0, 7 = 0) -» - (t/r) b , with — oc (-e) 7 . (23) 

T 

The entailed temporal power law, termed von Schweidler law, initiates the final decay 
of the correlators, which has a density and temperature independent shape <& q (i). In 
MCT, the (full) correlator thus takes the characteristic form of a two-step relaxation. 
The final decay, often termed a-relaxation, depends on e only via the time scale 
r(e) which rescales the time, t = t/r. Equation (22) establishes the crucial time 
scale separation between to and r, the divergence of r, and the stretching (non- 
exponentiality) of the final decay; it also gives the values of the exponents via A = 
r{l + b) 2 /r(l + 2b), and 7 = (l/a+l/6)/2. Using Eq. (17), the MCT-prediction for 
the divergence of the Newtonian viscosity follows [2, 37]. During the final decay the 
quiescent shear modulus also becomes a function of rescaled time, t = t/r, leading 
to rjo oc r(e); its initial value is given by the elastic constant at the transition, G^. 

The two asymptotic temporal power-laws of MCT also affect the frequency de- 
pendence of G" in the minimum region. The scaling function Q describes the mini- 
mum as crossover between two power laws in frequency. The approximation for the 
modulus around the minimum in the quiescent fluid becomes [37]: 



a + b 



+ 



(24) 



The parameters in this approximation follow from Eqs. (22,23) which give G m i n oc 
V^e and w m i n oc (— e) 1 / 2 ". Observation of this handy expression requires that the 
relaxation time r is (very) large, viz. that time scale separation holds (extremely well) 
for (very) small — e; even in the exemplary Fig. 9, the chosen distances to the glass 
transition are too large in order for the Eq. (24) to agree with the true /3-correlator Q, 
which is also included in Fig. 9. The reason for this difficulty is the aspect that the 
expansion in Eq. (20) is an expansion in \fs, which requires very small separation 
parameters s for corrections to be negligible. For packing fractions too far below the 
glass transition, the final relaxation process is not clearly separated from the high 
frequency relaxation. This holds in the experimental data shown in Fig. 6, where the 
final structural decay process only forms a shoulder. 
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On the glassy side of the transition, e > 0, the transient density fluctuations stays 
close to a plateau value for intermediate times which increases when going deeper 
into the glass, 



Entered into Eq. (17), the square-root dependence of the plateau value translates into 
the square-root anomaly of the elastic constant G^, and causes the increase of the 
yield stress close to the glass transition. 

Only, for vanishing shear rate, 7 = 0, an ideal glass state exists in the ITT ap- 
proach for steady shearing. All density correlators arrest at a long time limit, which 
from Eq. (25) close to the transition is given by <£ q (i — > 00, e > 0, 7 = 0) = f q = 
fq+h q y/e/{\ — X)+0(e). Consequently the modulus remains elastic at long times, 
g(t — > 00, e > 0,7 = 0) = Goo > 0. Any (infinitesimal) shear rate, however, melts 
the glass and causes a final decay of the transient correlators. The function Q initiates 
the decay around the critical plateau of the transient correlators and sets the common 
time scale for the final decay under shear 



Under shear all correlators decay from the plateau as function of \^t\ ; see e.g. Figs. 11, 
12, 21 and 22. Steady shearing thus prevents non-ergodic arrest and restores ergod- 
icity. Shear melts a glass and produces a unique steady state at long times. This con- 
clusion is restricted by the already discussed assumption to neglect aging of glassy 
states. It could remain because of non-ergodicity in the initial quiescent state, which 
needs to be shear-molten before ITT holds. Ergodicity of the sheared state, however, 
suggests aging to be unimportant under shear, and that it should be possible to melt 
initial non-ergodic contributions [8,9]. The experiments in model colloidal disper- 
sions reported in Sect. 3. 1 .2 and Sect. 6.2 support this notion, as history independent 
steady states could be achieved at all densities 5 . 

The described universal scenario of shear-molten glass and shear-thinnig fluid 
makes up the core of the MCT-ITT predictions derived from Eqs. (1 1) to (14). Their 
consequences for the nonlinear rheology will be discussed in more detail in the fol- 
lowing sections, while the MCT results for the linear viscoelasticity were reviewed 
in Sect. 3. Yet, the anisotropy of the equations has up to now prevented more com- 
plete solutions of the MCT-ITT equations of Sect. 2. Therefore, simplified MCT-ITT 
equations become important, which can be analysed in more detail and recover the 
central stability equations Eqs. (20,22). The two most important ones will be re- 
viewed next, before the theoretical picture is tested in comparison with experimental 
and simulations' data. 

5 An ultra-slow process causing the metastability of glassy states even without shear may 
have contributed to restore ergodicity in Refs. [31, 32]. 




(25) 




(26) 
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5 

Simplified models 

Two progressively more simplified models provide insights into the generic scenario 
of non-Newtonian flow, shear melting and solid yielding which emerge from the ITT 
approach. 



5.1 

Isotropically sheared hard sphere model 

On the fully microscopic level of description of a sheared colloidal suspension, affine 
motion of the particles with the solvent leads to anisotropic dynamics. Yet, recent 
simulation data of steady state structure factors indicate a rather isotropic distortion 
of the structure for Peo <C 1, even though the Weissenberg number Pe is already 
large [10,77]. Confocal microscopy data on concentrated solutions support this ob- 
servation [29]. The shift of the advected wavevector in Eq. (7) with time to higher 
values, intially is anisotropic, but becomes isotropic at longer times, when the magni- 
tude of q(t) increases along all directions. As the effective potentials felt by density 
fluctuations evolve with increasing wavevector, this leads to a decrease of friction 
functions, speed-up of structural rearrangements and shear-fluidization. Therefore, 
one may hope that an 'isotropically sheared hard spheres model' (ISHSM), which 
for 7 = exhibits the nonlinear coupling of density correlators with wavelength 
equal to the average particle distance (viz. the "cage-effect"), and which, for 7 7^ 0, 
incorporates shear-advection, captures some spatial aspects of shear driven decorre- 
lation. 



5.1.1 

Definition of the ISHSM 

Thus, in the ISHSM, the equation of motion for the density fluctuations at time t 
after starting the shear is approximated by the one of the quiescent system, namely 
Eq. (18a) (with r q = q 2 D s /S q ). The memory function also is taken as isotropic and 
modeled close to the unsheared situation [42, 76] 



(27a) 



with 



" Sn Sh Sri 



q • k c- k(t) + q p Cp (t ) [q • k c k + q • p c p 



(27b) 



where p = q k, and the length of the advected wavevectors is approximated by 
k(t) = fc(l + (t7/7 c ) 2 ) 1 / 2 and equivalently for p(t) . Note, that the memory function 
thus only depends on one time, and that shear advection leads to a dephasing of the 
two terms in the vertex Eq. (27b), which form a perfect square in the quiescent vertex 
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of Eq. (18c) without shear. This (presumably) also is the dominant effect of shear in 
the full microscopic MCT-ITT memory kernel in Eq. (14). The fudge factor 7 C is 
introduced in order to correct for the underestimate of the effect of shearing in the 
ISHSM 6 . 

The expression for the potential part of the transverse stress may be simplified to 

r°° u u t r ^ 5 c 'k S'i/ t \ 

cr = 7 / dtg(t,i) , with g(t,i) w dk — ^ *?. Jt) , (27c) 

1 Jo 60tt 2 J k{t)S 2 k fe(t) 

where, in the last equation Eq. (27c), the advected wavevector is chosen as k(t) — 
k(l + (^7) 2 /3) 1 / 2 , as follows from straight forward isotropic averaging of k(t) . For 
the numerical solution of the ISHSM for hard spheres using S q in PY approximation, 
the wavevector integrals were discretized as discussed in Sect. 3.1.1 and following 
Ref. [68], using M = 100 wavevectors from fc min = 0.1/R up to fc max = 19.9/i? 
with separation Ak = 0.2/ R. Again, time was discretized with initial step-width 
dt = 2 10~ 7 R 2 /D s , which was doubled each time after 400 steps [85]. The model's 
glass transition lies at (j) c — 0.51591, with exponent parameter A = 0.735 and w 
O.45/7 2 ; note that these values still change somewhat if the discretization is made 
finer. The separation parameter e = (<p — 4> c ) /4>c), and 7 are the two relevant control 
parameters determining the rheology. 

5.1.2 

Transient correlators 

The shapes of the transient density fluctuation functions can be studied with spatial 
resolution in the ISHSM. Figure 11 displays density correlators at two densities, 
just below (panel (a)) and just above (panels (b,c)) the transition, for varying shear 
rates. Panel (b) and (c) compare correlators at different wavevectors to exemplify the 
spatial variation. In almost all cases the shear rate is so small that the bare Peclet 
number Pe is negligibly small and the short-time dynamics is not affected. 

In the fluid case, the final or a-relaxation is also not affected for the two smallest 
dressed Peclet Pe values, but for larger Pe it becomes faster and less stretched; see 
the inset of fig. 11(a). 

Above the transition, the quiescent system forms an (idealized) glass [2, 37], 
whose density correlators arrest at the glass form factors f q from Fig. 10, and which 
exhibits a finite elastic constant G M , that describes the (zero-frequency) Hookian 
response of the amorphous solid to a small applied shear strain 7: a = Gooj for 
7^0; the plateau can be seen in Fig. 3 and, for intermediate times in Fig. 12. 
If steady flow is imposed on the system, however, the glass melts for any arbitrarily 
small shear rate. Particles are freed from their cages and diffusion perpendicular to 

6 Except for the introduction of the parameter 7 C , further quantitatively small, but qualita- 
tively irrelevant differences exist between the ISHSM defined here and used in Sect. 6.1 
according to Ref. [44], and the one originally defined in Ref. [42] and shown in Sect. 5.1; 
see Ref. [44] for a discussion. 
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Fig. 11. Normalized transient density correlators <£<j(t) of the ISHSM at wavevector q — 
3.4/d below (panel (a) at <j> — (f> c — 10 -4 ) and above (panel (b) at (f> = <p c + 10 -4 ) the 
transition for increasing shear rates Peo = 9 n * 10~ 8 with n = 0, . . . , 10 from right to left; 
the distances correspond to e = ±10 -3,53 . Curves for n = 9, 10 carry short and for n = 8 
long dashes; note the collapse of the two short dashed curves in (a). The insets show the data 
rescaled so as to coincide at <l>(t = r) = 0.1. Panel (c) shows glass correlators at another 
wavevector q — 7/d for parameters as labeled; from Ref. [76]. 
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Fig. 12. Transient non-Newtonian shear modulus g(t, 7) of the ISHSM in units of fcsT/c? 3 
for the packing fractions <f> — (f> c ± 1CP 3 (e = ±1CP 2,53 ; dashed/ dotted lines, respectively) 
for increasing shear rates Peo = 9" * 1CF 6 with n = 0, . . . , 8 from right to left; note the 
collapse of fluid lines for the smallest Pe ; from Ref. [76]. The solid line gives g(t, 7) for (f> = 
4> c and Peo = 1CP 6 . The inset shows the data rescaled so as to coincide at g(t = r, 7) = 5; 
note the collapse of the n — 6,7 & 8 curves for both (f> > 4> c and <f> < <f> c . 

the shear plane also becomes possible. Any finite shear rate, however small, sets a 
finite longest relaxation time, beyond which ergodicity is restored; see Figs. 1 l(b,c) 
and 12. 

The glassy curves at e > 0, panels (b,c), exhibit a shift of the final relaxation with 
Tj from Eq. (26) and asymptotically approach a scaling function (i/ T -y) • The mas- 
ter equation for the "yielding" scaling functions ^+ in the ISHSM can be obtained 
from eliminating the short-time dynamics in Eq. (18a). After a partial integration, 
the equation with d t 4> q (t) = is solved by the scaling functions: 

#+(t) = m+(t) I j* dP m+(t F) $+(F) , (28a) 
where t — t/r^, and the memory kernel is given by 

< w = ^ E v S (*) *t (*) *i-ki © • ( 28b > 

k 



While the vertex is evaluated at fixed shear rate, 7 = y -^f, it depends on the 
equilibrium parameters. The initialization for the correlator is given by 



*+(t^0) = /,, 



(28c) 
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with glass form factor taken from Eq. (21). The two-step relaxation and the shift of 
the final relaxation with r-y are quite apparent in Fig. 11. 

Figure 12 shows the transient shear modulus g(t,j) of the ISHSM from Eq. 
(27c), which determines the viscosity via r\ = J^°dtg(t 7 7). It is the time derivative 
of the shear stress growth function 77+ (i, -7) (or transient start up viscosity; here, 
the + labels the shear history [1, 87]), g(t, 7) = -^rj + (i, 7), and in the Newtonian- 
regime reduces to the time dependent shear modulus, g k (t). The g(t 7 7) shows all the 
features exhibited by the correlator of the density correlators in the ISHSM, and thus 
the discussion based on the stability analysis in Section 4 and the yielding scaling law 
carries over to it: g(t, 7) = G^, +h g Q + . ..; especially the dependence of the final 
relaxation step on rescaled time i/V-y is apparent in the glass curves. But in contrast 
to the density correlators & q (t), the function g(t,j) becomes negative (oscillatory) 
in the final approach towards zero, an effect more marked at high Pe. This behavior 
originates in the general expression for g(t, 7), Eq. (27c), where the vertex reduces 
to a positive function (complete square) only in the absence of shear advection. A 
overshoot and oscillatory approach of the start up viscosity to the steady state value, 
r] + (t — > 00, 7) — > 77(7), therefore are generic features predicted from our approach. 

5.1.3 

Flow curves 

As discussed in Section 4, in the fluid, MCT-ITT finds a linear or Newtonian regime 
in the limit 7^0, where it recovers the standard MCT approximation for New- 
tonian viscosity 770 of a viscoelastic fluid [2,37]. Hence a —* 7770 holds for Pe 
<C 1, as shown in Fig. 13, where Pe calculated with the structural relaxation time 
r is included. As discussed, the growth of r (asymptotically) dominates all trans- 
port coefficients of the colloidal suspension and causes an proportional increase in 
the viscosity 77. For Pe > 1, the non-linear viscosity shear thins, and a increases 
sublinearly with 7. The stress versus strain rate plot in Fig. 13 clearly exhibits a 
broad crossover between the linear Newtonian and a much weaker (asymptotically) 
7-independent variation of the stress. In the fluid, the flow curve takes a S-shape in 
double logarithmic representation, while in the glass it is bent upward only. 

Above the transition, the glass melts for any shear rate. Nonetheless, a finite 
limiting stress (yield stress) must be overcome in order to maintain the flow of the 
glass: 

<y{i, £ > 0) > <7 + (e > 0) = lim a(j, e > 0) . 

7— >o 

For e > and 7^0, the time for the final decay , Eq. (26), can become 
arbitrarily slow compared to the time characterizing the decay onto f q . Inserting the 
scaling functions <P + from Eq. (28) into the expression Eq. (27c) for the stress, the 
long time contribution separates out. Importantly, the integrands containing the 
functions depend on time only via t — t/r^ oc jt, so that nontrivial limits for the 
stationary stress follow in the limit 7 — > 0. In the ISHSM for e > 0, the yield stress 
is given by: 
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Fig. 13. Steady state shear stress a in units of ksT/d 3 versus Peo = jd 2 /Do, for the 
ISHSM at various distances from its glass transition, <f> — cj> c as labeled; circles correspond 
to fluid, diamonds to the critical, and squares to glassy densities; from Ref. [76], where the 
additional lines are discussed. For the fluid cases, cj> < (j> c , dashed lines indicate Newtonian 
fluid behavior, a — rjj, while vertical bars mark Pe= jr = 1, with the structural relaxation 
time taken from <P q = 7/d (t = t) = 0.1. The additional stress which would arise from the 
background solvent viscosity, a — "/r/oo, is marked by a dot-dashed line. For the critical 
density, <f> c , the critical yield stress, cr+ = 6.04, is shown by a horizontal bar. The inset shows 
the rise of the dynamical yield stress a + = a(e > 0, 7 — > 0+) in the glass together with a 
fitted power-law asymptote, a + = cr+ + 112V</> — <t>c- 



°+ = fit fdk Z% U+ it)) 2 , (28d) 

60tt 2 J J k{t)S 2 k V fe W V V 

The existence of a dynamic yield stress in the glass phase is thus seen to arise from 
the scaling law in Eq. (28), which is clearly borne out in the Figs. 11 and 12. The 
yield stress arises from those fluctuations which require the presence of shearing to 
prevent their arrest. Even though cr+ requires the solution of dynamical equations, 
in MCT-ITT it is completely determined by the equilibrium structure factor S q . This 
may suggest a connection of MCT to the potential energy paradigm for glasses, as 
recently discussed [78,79]. One might argue that a + arises because the external 
driving allows the system to overcome energy barriers so that different metastable 
states can be reached. This interpretation would agree with ideas from spin-glass [9] 
and soft-glassy rheology [6-8]. MCT-ITT indicates how shear achieves this in the 
case of colloidal suspensions. It pushes fluctuations to shorter wavelengths where 
smaller particle rearrangements cause their decorrelation. 

The increase of the amplitude of the yielding master functions <P + in Eq. (28) 
originates in the increase of the arrested structure in the unsheared glass Eq. (25). 
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In consequence the yield stress should rapidly increase as one moves further into 
the glass phase cr+ - cr+ oc \fe should (approximately) hold; see however [80] 
for the more complicated rigorous expression. Indeed, the inset of Fig. 12 shows a 
good fit of this anomalous increase to the numerical data. It is one of the hallmarks 
of the weakening of the glass upon approaching the glass transition from the low 
temperature or high density side. 



5.2 

Schematic F$ -model 

The universal aspects described in the previous Section 4 are contained in any ITT 
model that contains the central bifurcation scenario and recovers Eqs. (20,22). Equa- 
tion (20) states that spatial and temporal dependences decouple in the intermediate 
time window. Thus it is possible to investigate ITT models without proper spatial 
resolution. Because of the technical difficulty to evaluate the anisotropic functionals 
in Eqs. (1 ld,14), it is useful to restrict the description to few or to a single transient 
correlator. The best studied version of such a one-correlator model is the R, , -model. 



5.2.1 

Definition and parameters 

In the schematic -model [76] a single 'typical' density correlator &(t), conve- 
niently normalized according to <P(t — > 0) = 1 — rt, obeys a Zwanzig-Mori memory 
equation which is modeled according to Eq. (13) 

d t ${t) + r |*(t) + J dt' m(t - t') d t >${t') J = . (29a) 

The parameter r mimics the short time, microscopic dynamics, and depends on 
structural and hydrodynamic correlations. The memory function describes stress 
fluctuations which become more sluggish together with density fluctuations, because 
slow structural rearrangements dominate all quantities. A self consistent approxima- 
tion closing the equations of motion is made mimicking Eq. (14a). In the F$ -model 
one includes a linear term (absent in Eq. (14a)) in order to sweep out the full range 
of A values in Eq. (21), and in order to retain algebraic simplicity: 

l + (W7c) 2 

This model, for the quiescent case 7 = 0, was introduced by Gotze in 1984 
[37,81] and describes the development of slow structural relaxation upon increasing 
the coupling vertices Vi > 0; they mimic the dependence of the vertices in Eq. 
(14b) at 7 = on the equilibrium structure given by S q . Under shear an explicit 
time dependence of the couplings in m(t) captures the decorrelation by shear in Eq. 
(14b). The parameter j c sets a scale that is required in order for the accumulated 
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strain jt to matter. Shearing causes the dynamics to decay for long times, because 
fluctuations are advected to smaller wavelengths where small scale Brownian motion 
relaxes them. Equations (29a,29b) lead, with <P(t) = f c + (1 - f c ) 2 Q(t, e, 7), and the 
choice of the vertices v 2 = v% = 2, and v\ = + e (1 — f°) // c , where vf — 0.828, 
to the critical glass form factor f c = 0.293 and to the stability equation (22), with 
parameters 

A = 0.707 , c (7) = 0.586/ 7c 2 , and t = 0A26/T . 

The F$ -model possesses a line of glass transitions where the long time limit / = 
<P(t — ► 00) jumps discontinuously; it obeys the equivalent equation to Eq. (21). The 
glass transition line is parameterized by (v^v^) — ((2A — 1), 1)/A 2 with 0.5 < A < 
1, and f c = 1 — A. The present choice of transition point (vf, v%) is a typical one, 
which corresponds to the given typical A-value. The separation parameter e is the 
crucial control parameter as it takes the system through the transition. 

For simplicity, the quadratic dependence of the generalized shear modulus on 
density fluctuations is retained from the microscopic Eq. (1 Id). It simplifies because 
only one density mode is considered, and as, for simplicity, a dependence of the 
vertex (prefactor) v a on shear is neglected 

g(t,j) =v a $ 2 (t)+r ]oo 5{t-0+) . (29c) 

The parameter 7700 characterizes a short-time, high frequency viscosity and models 
viscous processes which require no structural relaxation, like in the general case Eq. 
(15). Together with r, it is the only model parameter affected by HI. Steady state 
shear stress under constant shearing, and viscosity then follow via integrating up the 
generalized modulus: 



r- oc r oc 

(7 = 777 = 7 / dtg(t)= A / / dt v a $ 2 (t) + 7 r\ 
Jo Jo 



(30) 



Also, when setting shear rate 7 = in Eqs. (29a,29b), so that the schematic correlator 
belongs to the quiescent, equilibrium system, the frequency dependent moduli are 
obtained from Fourier transforming: 

/>oo 

G'(uj)+iG"(uj) = iuj dte- tuJt v a $ 2 {t)\. +101^ . (31) 
Jo 7 

Because of the vanishing of the Fourier-integral in Eq. (31) for high frequencies, the 
parameter 7700 can be identified as high frequency viscosity: 

lim G"(u)/u = 77^ , with vl=Voo- (32) 

At high shear, on the other hand, Eq. (29b) leads to a vanishing of m(t), and Eq. (29a) 
gives an exponential decay of the transient correlator, <P(t) — > e~ rt for 7^0. The 
high shear viscosity thus becomes 

VL = . lim o-(i)/7 = Voc + ^ = C + ^ ■ (33) 
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5.2.2 

Correlators and stability analysis 

Representative solutions of the -model are summarized in Fig.9, which bring out 
the discussed universal aspects included in all ITT models. For small separation pa- 




Fig. 14. Numerically obtained transient correlators $ (t) (solid lines) for e = 0.01 (red, 1), 
e = (orange, 2), e = -0.005 (green, 3) and e = -0.01 (blue, 4) for the F$ -model from 
Ref. [80]. All curves were calculated with j — 10~ 7 . The dashed lines show the corresponding 
numerically obtained functions f c + (1 — f c ) 2 Q (t). 



rameters and shear rates the correlators develop a stretched dynamics located around 
the critical plateau value / c , according to Eqs. (20) to (22). The discussion of the dy- 
namics around this plateau was topic of Section 4. Figure 14 shows these aspects in 
the -model and presents typical correlators and the corresponding /3-correlators. 



Nonlinear rheology 



45 



The latter describe how the glassy structure, which is present on intermediate times, 
is molten either because the density is too low, or the temperature too high, or, al- 
ternatively, because of the effect of shearing. For long times Q (t) merges into the 
linear asymptote —t/r^ from Eq. (26). In the liquid region, for short times Q (t) fol- 
lows (t/to ) a from Eq. (22b) and merges into the second power law — (t /to ) b from 
Eq. (23) for intermediate times with the von MCT Schweidler exponent [37]. In the 
transition region close to e = 0, after following (t/to) ° the function Q (t) merges 
directly into the long time asymptote — t/r-y. In the yielding glass region, Q (t) fol- 
lows (t/to)~ a , arrests on the plateau value y/e/ (1 — A) for intermediate times and 
merges into the linear asymptote —t/r^ only for long times. So we can summarize 
that the short- and long time asymptotes are common for all e if 7 7^ is common. 
Fig. 15 shows an overview of the properties of Q (t). 

The present /3-scaling law bears some similarity to the one presented by Gotze 
and Sjogren for the description of thermally activated processes in glasses [54,82]. In 
both cases, ideal glass states are destroyed by additional decay mechanisms. Yet, the 
ITT equations and the generalised MCT equations differ qualitatively in the mecha- 
nism melting the glass. The similarity between both scaling laws thus underlines the 
universality of the glass stability analysis, which is determined by quite fundamental 
principles. In Eq. (22), the shear rate can only be a relevant perturbation (at long 
times) if it appears multiplied by time itself. Symmetry dictates the appearance of 
(jt) 2 , because the sign of the shear rate must not matter. The aspect that shear melts 
the glass determines the negative sign of (jt) 2 . 

The melting of the glassy structure during the yield process of Eq. (28) can be 
explicitly evaluated in the schematic F$ -model at e = 0. The yield master function 
does not depend on 7, and while its form is model-dependent, its initial decay follows 
from the universal stability equation Eq. (22). Figure 16 shows numerical results, 
which can be well approximated by an exponential function. 

The qualitative agreement between the transient correlators of the ISHSM and the 
schematic -model support the simplification to disregard the spatial structure. 
The only cost to be paid, is the fixed plateau value f c , which can be varied with 
wavevector in the ISHSM, but not in the schematic model. 

5.2.3 

Asymptotic laws of flow curves 

(7) 

A major advantage of the simplified F{ 2 ' -model is that it allows for asymptotic ex- 
pansions that qualitatively capture the flow curves. They can thus be investigated 
in detail addressing such questions as e.g. for the existence of power-law shear thin- 
ning [1], or the dependence of the yield stress on separation parameter. Fig. 17 shows 
an overview of the numerically obtained flow curves and the corresponding asymp- 
totic results given by the so-called /1-formula. While the glass flow curves exhibit 
an upward curvature only, the fluid curves show a characteristic S-shape, where the 
initial downward curvature changes to an upward one for increasing shear rate. Both 
behaviors are captured by the asymptotic expansions. For positive separation param- 
eters the range of validity of the /1-formula is given by |e| <C 1 and |-yto| <C 1. These 
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Fig. 1 5. An overview of the properties of Q (t) (solid lines) for the same values for e and 7 as 
in Fig. 14; from Ref. [80]. The dotted lines show the leading asymptotes for the corresponding 
time scales: The critical decay (t/to)" (a), the von Schweidler law — (i/Vo) 6 (b), the arrest 
on the plateau value y/e/ (1 — A) (c) and the shear-induced linear asymptote —t/rj (d). The 
dashed line shows a generalization of the latter law evaluated to higher order with a fitted 
parameter ai = -1.39 • 10 3 (at e = 0). 



two requirements ensure that Q (t) describes the dynamics of <P (t) with a sufficiently 
high accuracy; see Figs. 14 and 15. For sufficiently small negative separation param- 
eters the yl-formula is valid in finite shear rate windows only, as it does not reproduce 
the linear asymptotes for low shear rates. Precise criteria for the range of its valid- 
ity are known [80], and Fig. 17 presents an overview of the flow curves and their 
asymptotic laws. 

While the detailed discussion of the flow curves and their asymptotics leads be- 
yond the present review, see Ref. [80], the important conclusions from Fig. 17 in the 
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Fig. 1 6. The numerically determined master function of the yielding-process obeying a 'time- 
shear-superposition principle' (dotted line) $ + {t) = <Po (t) (label 0); from Ref. [80]. The 
solid lines show numerically obtained transient density correlators for e = and 7 = 1CP 7 
(red, a), 7 = 10" 9 ( green, b) and 7 = 10 12 (blue, c), plotted as functions of the rescaled 
time t. The plots demonstrate that the rescaled correlators converge to the yield master function 
So (t) from the analog of Eq. (28) in the -model for 7^0, the blue curve (c) is already 
quite close to the master curve (0). The dashed line shows <Po (t) + ai liT &i (t) for 7 = 
10~ r and the same numerical value for ai as in Fig. 15, using the leading correction 4>i (tj 
(1), which is shown in the lower panel. This first order expansion already describes quite well 
the shear-induced decay of the red curve (a). The inset demonstrates that the master function 
^0 (i) (dotted line) can be well approximated by an exponential function (solid line); the 
curves overlap completely. 
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Fig. 17. Overview of the numerically obtained flow curves (solid lines) and the asymptotic 
/1-formula evaluated numerically (dashed lines); from Ref. [80]. The liquid curves in the upper 
panel are shown for e = — 10~ 7 (red, 1), e = — 1CP 6 (orange, 2), e = — 1CP 5 (green, 3) and 
e = — (blue, 4). The lower panel shows the glassy curves for e = (red, 1), e = 10~ J 
(orange, 2), e = 10 -4 (green, 3) and e = 10 -3 (blue, 4). Crosses mark the points with 

2a 

\e\ = s-y = \jto \ i+o . The natural upper boundary for the shear rate, where the range 
of validity of the yl-formula is left, is also indicated. For e < 0, the natural lower limits for 
the shear rates, below which the /1-formula does not describe the flow curves, are marked by 
diamonds. 
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Fig. 18. The upper panel shows numerically obtained flow curves (solid lines) for e = 
-10~ 3 (red, 1), e = -1CT 2 (green, 2) and e = -1CT 1 (blue, 3). The dotted lines show the 
corresponding inflection tangents, with exponents p — 0.16, 0.35, and 0.63 from left to right. 
The dashed line shows the numerically evaluated yl-formula for e = — 10 -3 . The shear rate 
with e = — e-y is marked by a cross. The lower panel shows the corresponding results for the 
viscosity; fromRef. [80]. 



present context are: that the universal aspects discussed in Sect. 4 are recovered, that 
qualitative agreement is obtained with the results of the ISHSM, and that analytical 
expressions for the flow curves can be obtained. For example, the critical flow curve 
follows a generalized Herschel-Bulkley law: 

3 

a (e = 0, 7) = cr+ ^2 c„ I7/7* P™ , 

n=0 

where <j+ is the critical dynamic yield stress and 7* defines a natural scale for the 
shear rates in the asymptotic expansion; the upper limit 3 of the summation is dis- 
cussed in Ref. [80]. At the transition, this law describes the flow curve correctly 
for sufficiently small shear rates, see Fig. 17. This result also implies a generalized 
power-law weakening of the yield stress a + (e) when approaching the glass transition 
for e \ 0+, which is shown in the inset of Fig. 13. In fluid states for e < 0, the flow 
curves in double logarithmic presentation, viz. log 10 (a) as function of log 10 (7), 
show an inflection point defined by 

<f (log 10 (a)) _ Q 
d (login (7)) 2 

But then in some finite shear rate windows the flow curves can be approximated by 
the corresponding inflection tangents. The slopes p of the inflection tangents can be 
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interpreted as exponents occurring in some pseudo power laws: 

a oc 7 P . <^> 77 oc . 

Figure 18 shows some examples. The asymptotic formula also describes the neigh- 
borhood of the inflection point correctly for sufficiently small e, but does not rep- 
resent a real power law. In the framework of asymptotic expansions, there is thus 
no real exponent p. The power-law shear thinning often reported in the literature, in 
the ITT-flow curves thus is a trivial artifact of the double logarithmic plot. Rather, 
the flow curves on the fluid side exhibit a characteristic S-shape. While this shape 
is rather apparent when plotting stress versus shear rate, plotting the same data as 
viscosity as function of shear rate hides it, because the vertical axis gets appreciably 
stretched. 




Fig. 19. Reduced flow curves for a core-shell dispersion at an effective volume fraction of 
4> e ff = 0.580; data from Ref. [32], analysis from Ref. [80]. Here Rh denotes the hydrody- 
namic radius and Do the self diffusion coefficient of the colloidal particles; fcsT is the thermal 
energy. The solid line (red) shows the result for the fitted i 7 ^' -model with «| = 2.0. The fit- 
ted parameters are: e = -0.00042, -y c = 0.14, v„ = 70k B T/R 3 H , F = 80D /R 2 H and 
t?oo = 0.394fesT/i?jjDo. The dashed line shows the corresponding result for the yl-formula. 
The dotted line shows the inflection tangent of the numerically determined flow curve with a 
slope of p — 0.12. The inset shows the corresponding results for the viscosity. 



5.2.4 

Test of asymptotics in a polydisperse dispersion 

While the asymptotic expansions in the previous section provide an understanding 
of the contents of the MCT-ITT scenario, experimental tests of the asymptotic laws 
require flow curves over appreciable windows in shear rate. 
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Fig. 20. Reduced flow curves for a core-shell dispersion at an effective volume fraction of 
4> e ff = 0.629; quantities as defined in the caption of Fig. 19. The solid line (red) shows the 
result for the fitted F$ -model with «| = 2.0. The fitted parameters are: s = 0.000021, 
7 C = 0.16, v a = 115k B T/R 3 H , T = 120D /R 2 H and ^ = 0.431fc B T '/ R H D . The 
dashed line shows the corresponding result for the /1-formula. The dotted line shows the fitted 
Herschel-Bulkley law given by Eq. (34) with the analytically calculated exponent m = 0.489; 
data from Ref. [32], analysis from Ref. [80]. 

Figure 19 and Fig. 20 show experimental data recently obtained by Siebenburger 
et al. [32] on polydisperse dispersions of the thermosensitive core-shell particles in- 
troduced in Sect. 3.1.2 [30]. In all cases stationary states were achieved after shearing 
long enough, proving that ageing could be neglected even for glassy states. Because 
of the appreciable poyldispersity in particle size (standard deviation 17 %) crystal- 
lization could efficiently be prevented and flow curves over extremely wide windows 
could be obtained. Two flow curves from their work can be used to test the asymp- 
totic results. 

Fig. 19 shows the result for a liquid-like flow curve where the asymptotic A- 
formula holds for approximately four decades. The pseudo power law resulting from 
the inflection tangent of the flow curve holds for approximately two decades within 
the range of validity of the yl-formula. 

Fig. 20 shows the result for a flow curve, where a small positive separation pa- 
rameter was necessary to fit the flow curve and the linear viscoelastic moduli simul- 
taneously. The data are compatible with the (ideal) concept of a yield stress, but fall 
below the fit curves for very small shear rates. This indicates the existence of an 
additional decay mechanism neglected in the present approach [31,32]. Again, the 
yl-formula describes the experimental data correctly for approximately four decades. 
For higher shear rates, an effective Herschel-Bulkley law 



a{% < | 7 | < r,e 



0) = CTQ + <Tl \jtf)\ 



in 



(34) 
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with constant amplitudes and exponent m = 0.49 can be fitted in a window of ap- 
proximately two decades. The constant <r is not the actual yield stress, er+, which is 
obtained in the limit of vanishing shear rate, er + = <r(j — > 0), but is larger, 

The experimental data of the polydisperse samples, which exhibit structural dy- 
namics over large windows, and their fits with the full schematic model, will be taken 
up again in Sect. 6.2, where additionally the linear response moduli are considered, 
as had been done Sect. 3.1.2 for the less polydisperse samples affected by crystal- 
lization. 

6 

Comparison of theory and experiment 

As MCT-ITT contains uncontrolled approximations, justification to studying it, can 
be obtained only from its power to rationalize experimental observations. Because 
the transient density fluctuations are the central quantity in the approach, density 
correlators shall be considered first. Flow curves have been studied in most detail 
experimentally and in simulations, and thus are considered next. 

6.1 

ISHSM and single particle motion under steady shear 

Detailed measurements of the stationary dynamics under shear of a colloidal hard 
sphere glass have recently been obtained by confocal microscopy [29]. Single parti- 
cle motion was investigated in a shear-molten glass at roughly the wavevector inverse 
to the average particle separation. Figure 21 shows self-intermediate scattering func- 
tions measured for wavevectors along the vorticity direction where neither affine 
particle motion nor wavevector advection appears. The stationary correlators deep 
in the glass, for shear rates spanning almost two decades, are shown as function of 
accumulated strain jt, to test whether a simple scaling r-y ~ I/7 as predicted by 
Eq. (28) holds. Small but systematic deviations are apparent which have been inter- 
preted as a power law r-y ~ j~ - 8 [29,83]. ISHSM computations were performed 
for a nearby wavevector where S q is around unity so that coherent and incoher- 
ent correlators may be assumed to be similar [88]. Additionally, for the comparison 
it was assumed that time dependent transient and stationary fluctuation functions 
agree. The yielding master function from Eq. (28) in ISHSM can be fitted to the data 
measured at small effective Peclet numbers Pe e ff, by using for the phenomenologi- 
cal 'strain rescaling parameter' 7 C = 0.033; the smallness of the fitted value, which 
would be expected to be of order unity, is not yet understood. The effective Peclet 
number Pe cff = AR 2A f/D s with D s /D = 0.1 taken from [15] measures the impor- 
tance of shear relative to the Brownian diffusion time obtained from the short time 
self diffusion coefficient D s at the relevant volume fraction. At the larger effective 
Peclet numbers, Pe c ff > 0.5, for which the short-time and final (shear-induced) re- 
laxation processes move closer together, the model gives quite a good account of the 
7-dependence. 
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Fig. 21 . Steady state incoherent intermediate scattering functions & q (t) as functions of accu- 
mulated strain -yt for various shear rates 7; the data were obtained in a colloidal hard sphere 
dispersion at packing fraction (f> = 0.62 (at e m 0.07) using confocal microscopy [29]; the 
wavevector points in the vorticity (z) direction and has q = 3.8/ R (at the peak of S q ). The 
effective Peclet numbers Pe c g = 4R 2 -y/D s are estimated with the short time self diffusion 
coefficient D s m Do/ 10 at this concentration [15]. ISHSM calculations with separation pa- 
rameter e — 0.066 at qR = 3.9 (PY-S q peaking at qR = 3.5), and for strain parameter 
7 C = 0.033, are compared to the data for the Pe c ff values labeled. The yielding master func- 
tion at Pe c ff = lies in the data curves which span 0.055 < Pe c fr < 0.45, but discussion of 
the apparent systematic trend of the experimental data would require ISHSM to better approx- 
imate the shape of the final relaxation process; from Ref. [44]. 

The shape of the final relaxation step in a shear-molten glass can be studied even 
more closely in recent computer simulations, where a larger separation of short and 
long time dynamics could be achieved [84]. In these molecular dynamics simula- 
tions of an undercooled binary Lenard- Jones mixture, schematic ITT models give a 
good account of the steady state flow curves, 17(7) [77, 89]; this will be discussed in 
Sect. 6.3. Figure 22 shows the corresponding stationary self intermediate scattering 
functions for a wavevector near the peak in S q , oriented along the vorticity direction, 
for shear rates spanning more than four decades. Collapse onto a master function 
when plotted as function of accumulated strain is nicely observed as predicted by 
Eq. (28). At larger shear rates, the correlators raise above the master function; this 
resembles the behaviour observed in the confocal experiments in Fig. 21, and in the 
theoretical calculations in Figs. 11 and 16. Assuming again that transient coherent 
correlators can be fitted to stationary incoherent ones, the shape of the master func- 
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Fig. 22. Steady state incoherent intermediate scattering functions <Pq(t) measured in the vor- 
ticity direction as functions of accumulated strain jt for various shear rates 7; data from 
molecular dynamics simulations of a supercooled binary Lenard-Jones mixture below the 
glass transition are taken from Ref. [84]. These collapse onto a yield scaling function at long 
times. The wavevector is q = 3.55/i? (at the peak of S q ). The quiescent curve, shifted to 
agree with the one at the highest 7, shows ageing dynamics at longer times outside the plot- 
ted window. The apparent yielding master function from simulation is compared to the ones 
calculated in ISHSM for glassy states at or close to the transition (separation parameters e as 
labeled) and at nearby wave vectors (as labeled). ISHSM curves were chosen to match the 
plateau value f q , while strain parameters 7 C = 0.083 at e — (solid line) and 7 C = 0.116 at 
e = 10~ 3 (dashed line) were used; from Ref. [44]. 



tion can be fitted with the ISHSM, using again an unaccountedly small strain param- 
eter 7 C . After this rescaling, modest but visible differences in the shapes remain: the 
theoretical master function decays more steeply than the one from simulations. 

Overall, theory and experiment agree in finding a two step relaxation process, 
where shear has a strong effect on the final structural relaxation, while the short 
time diffusion is not much affected. This supports the central MCT-ITT prediction 
that shearing speeds up the structural rearrangements in a concentrated dispersion 
close to vitrification. More detailed comparisons await better theoretical calculations 
where the effect of shear on the stationary density fluctuation functions is taken into 
account more faithfully than in the ISHSM. 
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6.2 

F^-model and shear stresses in equilibrium and under flow in a 
polydisperse dispersion 

A central result of MCT-ITT concerns the close connection between structural re- 
laxation at the glass transition and the rheological properties far from equilibrium. 
The ITT approach aims to unify the understanding of these two phenomena, which 
were introduced wrt. experimental data in Sects. 3.1.1 and 5.2.4, respectively. MCT- 
ITT requires, as sole input, information on the equilibrium structure (namely S q ), 
and, first gives a formally exact generalization of the shear modulus to finite shear 
rates, g(t,j), which is then approximated in a consistent way. A novel dense col- 
loidal dispersion serves as experimental model system, whose linear and nonlinear 
rheology can be determined over very broad windows of control parameters. The 
generalized modulus g(t,j) can thus be investigated as function of shear rate and 
time (more precisely frequency), and the MCT-ITT approach can be tested excruci- 
atingly. Thermosensitive core-shell particles consisting of a polystyrene core and 
a crosslinked poly(N-isopropylacrylamide)(PNIPAM) shell were synthesized and 
their slightly polydisperse dispersions (standard deviation 17 %) characterized in 
detail [32]; see Sect. 3.1.1. While their precise structure factor has not been mea- 
sured yet, the system can be well considered a slightly polydisperse mixture of hard 
spheres. Because polydispersity prevents crystallization and the lack of attractions 
prevents demixing and coagulation, this system opens a window on structural re- 
laxation, which can nicely be tuned by changing the effective packing fraction by 
varying temperature. 

Shear stresses measured in non-linear response of the dispersion under strong 
steady shearing, and frequency dependent shear moduli arising from thermal shear 
stress fluctuations in the quiescent dispersion were measured and fitted with results 
from the schematic F-j^ -model. Some results from the microscopic MCT for the 
equilibrium moduli were included also; see Fig. 23. The fits with quiescent MCT 
for (monodisperse) hard spheres using the PY S q support the finding of Sect. 3.1.2 
that MCT accounts for the magnitude of the stresses at the glass transition semi- 
quantitatively. Because of the polydispersity of the samples, which is neglected in 
the calculations performed according to the presentation in Sect. 3.1.2, somewhat 
larger rescaling factors c y are required; they are included in Table 1 . Also, the critical 
packing fraction <j) c of the glass transition again is reproduced with some small error. 
Because of polydispersity, the experimental estimate </£ ff w 0.625 [32] lies somewhat 
higher than in the (more) monodisperse case [31], which is as expected [86]. 

Figure 23 gives the comparison of the experimental flow curves and the linear 
response moduli G' and G" with theory for five given different effective volume 
fraction </> e ff adjusted by varying temperature. On the left-hand side the flow curves 
(7(7) are presented as function of the bare Peclet number Pe = kBT/(Qnr] s R^ { ) 7, 
on the right-hand side G' and G" are displayed as the function of the frequency- 
Peclet or Deborrah number Pe w = fc^T/ (GnrjsR^) lj, calculated with the respec- 
tive frequency u. Table 1 gathers the effective volume fractions together with the fit 
parameters of the F? 2 -model. Note that G' and G" have been obtained over nearly 7 
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orders of magnitude in frequency, while the flow curves extend over more than eight 
decades in shear rate. 

The generalized shear modulus g(t, 7) of the F± 2 -model (cf. eq. 29c) presents the 
central theoretical quantity used in these fits. Within the schematic model, the vertex 
prefactor v a is kept as a shear-independent quantity. It can easily be obtained from 
the stress and modulus magnitudes. Hydrodynamic interactions enter through 7700, 
which can be obtained from measurements done at high frequencies, and through 
r, which can be obtained via Eq. (33) from measurements done at high shear rates. 
Given these three parameters, both the shapes of the flow curves as well as the shapes 
of the moduli G' and G" may be obtained as function of the two parameters e and 
7/7 c . The former sets the separation to the glass transition and thus (especially) the 
longest relaxation time, while the latter tunes the effect of the shear flow on the flow 
curve. 

All measured quantities, namely a, G', and G" were converted to the respective 
dimensionless quantities by multiplication with R^/ksT where Rh is the hydro- 
dynamic radius at the respective temperature. As already discussed above, the ex- 
perimental control parameters 7 and u also were converted by Q-K^sR^j/ksT to the 
respective Peclet numbers. Evidently, both the reduced moduli, the Pe number, and 



0eff : v a r 7c 7700 £ mcro D s /D c y p 

l] r 1 



W-mo d. [%g] [jjj fe] micro 



0.530 -0.072000 18 20 0.0845 0.2250 -0.10 0.3 2.3 0.631 

0.595 -0.003500 48 50 0.1195 0.2400 -0.008 0.3 2.3 0.248 

0.616 -0.000420 70 80 0.1414 0.3938 -0.001 0.3 2.3 0.117 

0.625 -0.000170 85 90 0.1491 0.4250 -0.001 0.3 3.0 0.0852 

0.627 0.000021 115 120 0.1622 0.4313 0.002 0.3 3.5 - 

Table 1 . Packing fraction e ff and parameters v a , r, 7 C and r/oo of the fit using the schematic 
Fj 7 2 -model for the measurements shown in Fig. 23; from Ref. [32]. The parameters e, short 
time diffusion coefficient D a /Do, and rescaling factor c y from the microscopic linear re- 
sponse calculation using MCT , and the pseudo-power law exponent p are included, also. 



the packing fraction depend on the effective particle volume R? H . In the poly disperse 
sample, actually a distribution of values R? H exists, whose variance may be deter- 
mined by disc centrifugation at low concentration, and which is fixed for one given 
sample. Close to the glass transition the size distribution thus is (almost) density in- 
dependent, and R%(T) is the single experimental control parameter, whose small 
change upon varying temperature T drives the system through the glass transition. 

Figure 23 demonstrates that the rheological behavior of a non-crystallizing col- 
loidal dispersion can be modeled in a highly satisfactory manner by five parameters 
that display only a weak dependence on the effective volume fraction of the particles. 
Increasing the effective packing fraction drives the system towards the glass transi- 
tion, viz. e increases with </> e ff- Stress magnitudes (measured by v a ) also increase with 
e ff, as do high frequency and high shear viscosities; their difference determines r. 
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Fig. 23. Left column: Reduced flow curves (filled squares) for different 
volume fractions. The solid lines are the results of the schematic model, 
the dashed line represent the pseudo power law behaviour; from Ref. [32]. 
Right column: Reduced frequency dependent moduli for different volume fractions. 
Full symbols/solid lines represent G\ hollow symbols/dashed lines represent G". Thick lines 
are the results of the schematic model, the thin lines the results of the microscopic MCT 
Graphs in one row represent the continuous and dynamic measurements at one volume 
fraction, a and b at d> e « = 0.530 , c and d at d> e « = 0.595 , e and/ at 4>es = 0.616 , g and h at 



58 



Matthias Fuchs 



The strain scale 7 C remains around the reasonable value 10%. In spite of the smooth 
and small changes of the model parameters, the -F-^-model achieves to capture the 
qualitative change of the linear and non-linear rheology. The measured Newtonian 
viscosity increases by a factor around 10 5 . The elastic modulus G' at low frequen- 
cies is utterly negligible at low densities, while it takes a rather constant value around 
lOksT/R^ at high densities. An analogous observation holds for the steady state 
shear stress (7(7), which at high densities takes values around 0.3fcsT '/ 'R% when 
measured at lowest shear rates. For lower densities, shear rates larger by a factor 
around 10 7 would be required to obtain such high stress values. 

At volume fractions around 0.5 the suspension is Newtonian at small Peo- Ap- 
proaching the glass transition leads to a characteristic S-shape of the flow curves 
and the Newtonian region becomes more and more restricted to the region of small- 
est Peo- Concomitantly, a pronounced minimum in G" starts to develop, separating 
the slow structural relaxation process from faster, rather density independent mo- 
tions, while G' exhibits a more and more pronounced plateau. At the highest density 
(Fig. 23 panels i and j), the theory would conclude that a yielding glass is formed, 
which exhibits a finite elastic shear modulus (elastic constant) G^ = G'(u) — > 0), 
and a finite dynamic yield stress, a + = a(j — > 0). The experiment shows, however, 
that small deviations from this glass like response exist at very small frequencies and 
strain rates. Description of this ultra-slow process requires extensions of the present 
MCT-ITT which are discussed in Ref. [31]. 

Considering the low frequency spectra in G'(u) and G"(ui), microscopic MCT 
and schematic model provide completely equivalent descriptions of the measured 
data. Differences in the fits in Fig. 23 for Pe w < 1 only remain because of slightly 
different choices of the fit parameters which were not tuned to be close. These dif- 
ferences serve to provide some estimate of uncertainties in the fitting procedures. 
Main conclusion of the comparisons is the agreement of the moduli from micro- 
scopic MCT, schematic ITT model, and from the measurements. This observation 
strongly supports the universality of the glass transition scenario which is a central 
line of reasoning in the ITT approach to the non-linear rheology. 

6.3 

F^-model and flow curves of a simulated supercooled binary liquid 

In large scale molecular dynamics simulations a 80:20 binary mixture of Lennard- 
Jones (LJ) particles at constant density was supercooled under shear. This model has 
well known equilibrium properties and many aspects that can be understood con- 
sistenly within MCT [90]. To account for shearing, it was used together with Lees- 
Edwards boundary conditions and the SLLOD equations of motion to develop a lin- 
ear velocity profile. Note that in the simulation solvent effects are obviously lacking, 
and the simulated flow curves thus provide support for the notion that shear thinning 
can arise from shear-induced speed up of the structural relaxation; it is evident in 
Fig. 22. Because the microscopic motion is Newtonian, the theoretical description of 
this model goes beyond the framework of Sect. 2. Yet the universality of the struc- 
tural long time dynamics, predicted by MCT [70] and also retained in MCT-ITT, 
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Fig. 24. Flow curves 17(7) reaching from the supercooled to the glassy state of a 
simulated binary LJ mixture. The data points correspond to the temperatures T = 

0. 525, 0.5, 0.45, 0.44, 0.43, 0.42, 0.4, 0.38, 0.3, 0.2 and 0.01 in LJ-units (from bottom to 
top). .F-lj -model curves fitted by eye are included as lines. The inset shows the relation be- 
tween the fitted separation parameters and temperature. Units are converted by a = 1.5<r t heo 
and 7 = 1.37theo^; from [89]. The arrows mark the values of the extrapolated dynamic yield 
stresses a + (e) . 

supports to apply MCT-ITT to the simulation data. Moreover, the independence of 
the glassy dynamics on the employed microscopic motion was explicitly confirmed 
in simulations of the mixture [91]. This supercooled simple liquid has been char- 
acterized quite extensively under shear [10,49,77,84,89], and thus fits of the flow 
curves provide challenging tests to the schematic -model. 

Figure 24 shows the stress-shear rate dependence as flow curves, ranging from 
supercooled states to the glassy regime; LJ units are used as described in [77]. The 
solid lines are fits to the simulation data with the _F 12 -model, which reproduce 
the transition from a shear-thinning fluid to a yielding glass quite well. Coming 
from high shear rates, the flow curves of the supercooled state pass to the linear 
response regime in the lower left corner, indicated by a dashed-dotted line with slope 

1. On approaching the transition point, the linear response regime shifts to lower 
and lower shear rates. Beyond the fluid domain, the existence of a dynamic yield 
stress <j + = linx^o ct > is supported by the simulation results, which sustain 
a stress plateau over three decades in shear rate. Best -model fits are obtained 
for a T c = 0.4 suggesting a slightly lower transition temperature [92] as determined 
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Fig. 25. Dynamic yield stress estimated from the simulations of a supercooled binary LJ 
mixture under steady shear shown in Fig. 24, and its temperature dependence (in LJ units); 
from Ref. [77]. The estimate uses the stress values for the two lowest simulated shear rates, 
namely 7 = 10~ J (triangle) and 7 = 3 x 10 -6 (circle); the extrapolation with the i^j - 
model is shown by diamonds. At temperatures below T = 0.38, (almost) the same shear 
stress is obtained for both values of 7 and the extrapolation, indicating the presence of a yield 
stress plateau. 

from the simulations of the quiescent system, where T c = 0.435 was found [90]. 
The reason may be the ergodicity restoring processes which were also observed in 
the colloidal experiments shown in Figs. 7, 20, and 23 

The stress plateau is best developed for temperatures deep in the glassy phase 
extending over about two decades in shear rate. Its onset is shifted toward progres- 
sively lower 7 as the temperature is increased toward T c . This makes an estimate of 
the dynamic yield stress, cr + (T) = a(T; 7 — * 0), a difficult task for temperatures 
below but close to T c . Nevertheless, an estimate of a + (T) is interesting because it 
highlights the anomalous weakening of the glass when heating to T c . Testing the 
MCT predictions below T c has previously not been possible in simulations because 
of problems to reach the equilibrated or steady state at sufficiently low shear rates. 
Figure 25 significantly supports the notion of a glass transition under shear as it 
presents the first simulations result exibiting the predicted anomalous softening in 
an elastic property of the glass upon approaching the transition from the glass side, 
viz. upon heating. 

7 

Summary and outlook 

The present review explored the connection between the physics of the glass tran- 
sition and the rheology of dense colloidal dispersions, including in strong steady 
shear flow. A microscopic theoretical approach for the shear-thinning of concen- 
trated suspensions and the yielding of colloidal glasses was presented, which builds 
on the mode coupling theory (MCT) of idealized glass transitions. The extension 
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to strongly driven stationary states uses the so-called ITT (integration through tran- 
sients) approach, which leads to a scenario of shear melting a glass, whose universal 
aspects can be captured in simplified schematic models. Consecutive generalizations 
of ITT to arbitrary time-dependent states far from equilibrium [58] and to arbitrary 
flow geometries [93], have yielded a non-Newtonian constitutive equation applica- 
ble to concentrated dispersions in arbitrary homogeneous flows (not reviewed here), 
albeit still under the approximation that hydrodynamic interactions are neglected. 
Within the theory, this approximation becomes valid close to the glass transition and 
for weak but nonlinear flows, where the slow structural relaxation dominates the sys- 
tem properties, and where hydrodynamic interactions only affect the over-all time 
scale. 

The structural dynamics under flow is predicted to result from a competition be- 
tween local particle hindrance (termed cage effect) and the compression / stretching 
(i.e. advection) of the wavelength of fluctuations induced by the affine particle mo- 
tion with the flow. Measurements of the single particle motion in the stationary state 
under shear support the theoretical picture that shear speeds up the structural dy- 
namics, while instantaneous structural correlations remain rather unaffected. Model 
dispersions made of thermo-sensitive core-shell particles allow to investigate the 
close vicinity of the transition. Measurements of the equilibrium stress fluctuations, 
viz. linear storage and loss moduli, and measurements of flow curves, viz. nonlinear 
steady state shear stress versus shear rate, for identical external control parameters 
verify that the glassy structural relaxation can be driven by shearing and in turn itself 
dominates the low shear or low frequency rheology. 

In the employed theoretical approach, the equilibrium structure factor S q cap- 
tures the particle interactions. Theory misses an ultra-slow decay of all glassy states, 
and neglects (possible) ageing effects. 
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